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Preface 


This  thesis  examines  the  feasibility  of  placing  a  DIAL  lidar  on 
the  space  shuttle  ana  have  this  instrument  measure  the  intensity  of 
thunderstorms.  I  have  always  had  a  great  deal  of  interest  in 
meteorology  and  in  particular  the  thunderstorm.  As  an  undergraduate 
my  classmates  and  I  use  to  chase  thunderstorms  in  hopes  of  watching 
tornado  formation.  Now,  with  extensive  physics  and  remote  sensing 
courses  taught-  here  at  AFIT,  I  can  combine  my  meteorological  exper¬ 
ience  with  sensor  design  concepts  and  show  what  a  great  asset  the 
lidar  is  to  the  atmospheric  sciences. 

I  would  like  to  thank  my  advisor,  Major  Jim  Lange  (PhD) ,  who 
seems  to  always  be  two  steps  ahead  of  the  world,  for  all  of  his 
advice  and  support.  I  thank  the  AFIT  research  librarian,  Linda 
Stoddart,  for  locating  all  my  meteorological  references  which  re¬ 
quired  inter-library  loans.  Finally,  I  would  like  to  thank  my  wife. 
Gay,  for  eighteen  months  of  loving  support  during  very  trying  times 
here  at  AFIT. 


Craig  Z.  Lowery 


Abstract 


This  thesis  examines  the  feasibility  of  operating  a  lidar  (laser 
radar)  from  the  space  shuttle  and  having  this  instrument  measure  the 
intensity  of  thunderstorms.  Thunderstorm  intensity  is  monitored  by 
measuring  the  time  rate  of  change  in  temperature  of  the  top  0.5  km  of 
a  thunderstorm.  Severe  weather  occurring  on  the  ground  takes  place 
during,  or  just  after  a  period  of  rapid  cloud  top  cooling.  Tempera¬ 
ture  is  measured  with  two  wavelengths  using  the  differential  absorp¬ 
tion  lidar  (DIAL)  technique  which  determines  the  resonant  absorption 
of  oxygen  in  the  oxygen  "A"  band  near  .7700  pm.  One  wavelength  is 
set  at  .7696  pm  which  is  a  temperature  sensitive  oxygen  absorbing 
line  while  the  second  wavelength  is  set  at  .7614  pm  where  oxygen  ab¬ 
sorption  is  negligible.  The  lidar  can  measure  the  oxygen  resonant 
absorption  coefficient  at  the  heights  of  typical  thunderstorm  tops. 

A  temperature  value  is  recovered  by  placing  the  absorption  coefficient 
into  a  quickly  converging  iterative  expression.  Hypothetical 
thunderstorms  with  heights  from  10-17  km  are  probed  at  shuttle 
altitudes  ranging  from  100-250  km.  Success  of  the  system  is  based 
on  useful  values  of  the  signal-to-noise  ratio. 
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The  Space  Shuttle  Columbia  completed  its  fourth  test  flight  and 
became  operational  on  4  July  1982.  Now,  for  the  first  time,  a  space 
vehicle  can  carry  cargo  into  orbit  and  return  all  of  it  safely  back  to 
earth  in  a  relatively  economic  manner.  Such  a  magnificent  capability 
poses  an  interesting  question:  what  can  the  space  shuttle  be  used  for? 
Presentations  at  the  recent  Space  Symposium  held  at  the  U.S.  Air  Force 
Academy  on  11  October  1982  indicate  that  the  shuttle  is  an  excellent 
vehicle  to  conduct  scientific  research  from.  In  addition,  the  symposium 
revealed  that  proposed  scientific  experiments  have  yet  to  fully  utilize 
the  shuttle’s  unique  capabilities.  This  thesis  proposes  a  new  use  for 
the  space  shuttle. 

This  thesis  examines  the  feasibility  of  mounting  a  lidar  (laser 
radar)  in  the  shuttle's  cargo  bay  and  having  this  device  measure  the 
intensity  of  thunderstorms.  Thunderstorm  intensity  will  be  monitored 
by  measuring  the  time  rate  of  change  in  temperature  of  the  thunderstorm 
top.  Temperature  will  be  measured  using  a  recently  developed  DIAL  tech¬ 
nique  which  determines  the  resonant  absorption  of  oxygen  in  the  oxygen 
"A"  band  near  770  nm.  This  is  the  first  proposed  use  of  a  temperature 
sensing  laser  for  the  monitoring  of  thunderstorm  intensity. 

If  a  thunderstorm  becomes  severe,  a  great  deal  of  property  damage 
and  a  loss  of  life  is  possible.  Recent  research  indicates  that  severe 
weather  occurring  on  the  ground  took  place  during,  or  just  after  a  period 


of  rapid  cloud  top  cooling.  A  lidar  operating  from  the  space  shuttle 


could  measure  the  temperature  of  a  particular  thunderstorm  and  provide  an 


indication  of  the  storm's  intensity  or  severity.  A  shuttle- borne  lidar 
would  be  an  interesting  experiment  with  operational  possibilities  later 
on  such  as  providing  timely  warnings.  In  this  way,  such  a  sensor  could 
prevent  millions  of  dollars  of  damage  and  possibly  save  lives. 


II.  Literature  Review 


Characteristics  of  a  Thunderstorm 

A  thunderstorm  is  produced  by  a  cumulonimbus  cloud  and  is  always 
accompanied  by  thunder,  lightning,  rain,  wind,  and  occasionally  with 
hail  and/or  tornado  (s).  The  thunderstorm  and  its  associated  weather 
phenomenon  is  the  most  violent  storm  produced  in  the  earth's  atmosphere. 

In  the  continental  mid-latitudes,  thunderstorms  usually  occur  in  the 
spring  and  summer.  Those  storms  that  form  in  the  spring  are  typically 
triggered  by  synoptic  scale  (storms  1000-2500  km  across)  disturbances 
embedded  in  the  jet  stream.  These  spring-time  thunderstorms  can  be  quite 
intense.  Such  intense  or  severe  storms  can  reach  heights  greater  than 
14  km  and  produce  hail  (2  3cm  in  diameter),  strong  straight-line  winds 
( Jr  30  m  sec  ^),  intense  updrafts  (98  m  sec  ^  reported  by  Brown  and  Knupp, 
1980),  and  sometimes  tornados.  The  typical  summer-time  thunderstorm  is 
usually  less  severe  and  is  triggered  this  time  by  mesoscale  (storms  1- 
999  km  across)  disturbances  and/or  surface  heating. 

All  thunderstorms  can  occur  as  single  cells  by  themselves,  groups 
of  several  cells  producing  a  multicell  complex,  or  lined  up  singly  to 
form  a  squall  line.  Occasionally  a  single  local  storm  will  become  very 
severe  and  becomes  what  Browning  (1977)  refers  to  as  a  supercell.  A  few 
over-simplified  characteristics  of  severe  and  non-severe  thunderstorms 
can  be  found  in  Figure  1. 

The  shear  referred  to  in  Figure  1  is  the  horizontal  force  felt  by 
the  storm  as  wind  velocity  increases  with  height. 
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Figure  1.  Typical  thunderstorm  characteristics  (after  Rust  et  al,  1981 
Kropfli  and  Miller,  1976  with  chances).  " 


Internal  Structure  of  a  Thunderstorm 

Research  in  mesoscale  meteorology  is  dominated  by  investigations  of 
supercell  thunderstorms.  Such  research  is  carried  out  by  surface  observing 
networks,  radar,  aircraft,  and  satellite.  There  has  been  extensive  effort 
to  develop  computer  programs  which  simulate  thunderstorm  development; 
however,  to  this  day  the  internal  dynamics  of  a  thunderstorm  are  not  well 
understood,  and  with  the  lack  of  good  data,  computer  simulations  have  not 
verified  well  against  real  world  situations.  Thus,  most  thunderstorm 
models  are  of  a  conceptual  format  (Newton,  1963;  Browning,  1964;  Marwitz, 
1972  a.b.;  Lemon  and  Doswell,  1979b).  To  briefly  review  the  internal 
structure  of  a  severe  thunderstorm,  the  model  of  Lemon  and  Doswell  (1979b) 
is  chosen  since  it  represents  an  excellent  composite  of  current  research 
on  supercells. 

Figure  2  is  a  one-dimensional  view  of  a  supercell  storm  looking  down 
from  space.  The  heavy  dark  outline  represents  the  shape  of  the  storm  as 
viewed  from  a  weather  radarscope.  The  storm  is  moving  towards  the  north¬ 
east.  A  mesoscale  frontal  structure  is  formed  by  these  storms  which 
closely  resemble  synoptic  scale  cold  fronts  which  cover  entire  continents 
(Brooks,  1949). 

Figure  3  a-d  is  a  three-dimensional  schematic  of  a  supercell  with 
the  vertical  scale  distorted  and  features  above  9  km  omitted.  Figure  3a 
shows  the  main  updraft  rotating  upward.  The  forward  flank  downdraft  is 
formed  by  falling  precipitation  and  supplemented  by  mid-altitude  prevail¬ 
ing  winds  deflected  downward. 

Figure  3b  illustrates  the  unique  feature  of  a  supercell,  a  second 
downdraft.  This  rear  flank  downdraft  forms  as  cool  dry  air  slams  into  the 


Lemon  and  Doswell,  1979b) 


updraft  column  and  is  deflected  downward.  Figure  3c  shows  the  proposed 
location  of  a  tornado  existing  in  a  region  of  intense  vertical  velocity 
and  temperature  gradients. 

Figure  3d  shows  the  dissipating  stage  where  downdrafts  have  choked 
off  the  warm  moist  updraft.  Tornado  touchdown  is  possible  here. 

Top  Structure  of  a  Thunderstorm 

The  study  of  overshooting  cloud  tops  or  towers  protruding  above  the 
anvil  floor  by  several  kilometers  has  been  the  subject  of  much  research 
in  recent  years  (McCann,  1979).  Overshooting  tops  have  been  directly 
correlated  with  thunderstorms  producing  tornados,  hail,  and  strong  wind 
gusts  (see  Figure  1) . 

There  are  two  evolving  schools  of  thought  concerning  overshooting 
storm  tops.  The  first  theory  claims  that  severe  weather  will  occur  on 
the  ground  during  or  soon  after  an  overshooting  top  collapses.  Hard 
data  backs  this  theory.  The  second  theory  concerning  overshooting  tops 
claims  that  severe  weather  occurs  on  the  ground  during  or  soon  after 
rapid  top  growth.  Hard  data  backs  this  theory  also  (McCann,  1979) . 

Fujita  (1973a)  maintains  that  a  collapsing  top  will  create  a  tornado 
at  ground  level.  Here  it  is  claimed  that  the  updraft  carries  warm  moist 
air  above  the  anvil  floor.  This  now  protruding  top  soon  becomes  precipi¬ 
tation  heavy  causing  the  top  to  collapse  upon  itself.  The  resulting 
downdraft  twists  around  the  rotating  updraft  which  produces  enough  surface 
vorticity  to  generate  a  tornado. 

Lemon  at  al_  (1975)  claims  that  rotation  of  the  updraft  causes  the  top 
to  collapse  through  a  process  in  fluid  mechanics  called  vortex  valve.  Thi 


too  is  a  trigger  mechanism  for  tornadogcnesis. 


Adler  and  Fenn  (1979c)  studied  thunderstorm  vertical  growth  rates 
and  top  structure  from  infrared  geosynchronous  satellite  data.  These 
researchers  showed  that  the  rate  of  storm  growth  and  minimum  cloud  top 
temperature  correlate  with  reports  of  severe  weather  on  the  ground. 

Reynolds  (1980)  studied  top  temperatures  of  damaging  hailstorms  with 
infrared  geosynchronous  satellite  data.  He  determined  that  hailfall 
occurred  at  close  to  the  time  of  the  maximum  rate  of  storm  growth.  In 
other  words,  hail  was  reported  under  an  overshooting  top  as  the  top 
reached  its  maximum  height  (coldest  temperature). 


Forecasting  Thunderstorms 

The  National  Weather  Service  (NWS)  defines  a  severe  thunderstorm  as 
one  or  more  of  the  following:  3/4  inch  (1.9  cm)  or  larger  hail,  strong 
winds  of  50  knots  (2G  msec  ^)  or  greater,  and  tornado (s). 

Ostby  and  Higginbotham  (1982)  examined  tornado  predictability  and 
intensity.  They  considered  all  reported  tornados  (8,825)  for  the  period 
1971-1980.  Figure  4  shows  the  10  year  statistics  for  the  number  of  weak, 
strong,  and  violent  tornados  for  each  state.  Figure  5  shows  the  number 
of  tornado  deaths  during  the  10  year  period  which  occurred  in  either  a 


tornado  watch  or  in  a  severe  thunderstorm  watch. 


Lemon  (1979a)  points  out  that  the  NWS's  thunderstorm  and  tornado 
warning  programs  need  considerable  improvement.  The  Doppler  radar  will 
improve  warnings,  but  the  nation-wide  network  will  not  be  operational  until 
the  1990  time  frame  (Lemon  et^  al,  1977;  Burgess  and  Devore,  1979; 

Wilson,  1980).  Stereographic  observations  from  geosynchronous  satellites 
of  thunderstorm  height  changes  is  a  promising  new  tool;  however,  it  is 


only  now  being  demonstrated,  and  researchers  would  like  more  satellites 
(Hasler,  1981). 

Any  improvement  in  forecast  reliability  in  the  near  future  will 
lower  the  false  alarm  rate  (no  occurrence  of  a  tornado  when  forecast). 

A  high  false  alarm  rate  may  produce  public  apathy.  If  the  loss  of  life 
is  to  be  prevented,  severe  thunderstorm  and  tornado  warnings  must  improve 
(Johnson,  1977). 

Satellite  Observations  of  Thunderstorms 

Distinct  severe  thunderstorm  characteristics  were  first  identified 
by  using  visual  images  from  the  old  TIROS  1  satellite,  Whitney  and  Fritz 
(1961)  and  Whitney  (1963)  identified  severe  thunderstorms  to  be  (1)  vis¬ 
ibly  brighter  than  surrounding  clouds,  (2)  rounded  in  appearance,  (3) 
have  scalloped  edges,  (4)  some  evidence  of  cirrus  blowoff,  and  (5)  di- 
mensins  of  100  to  200  miles.  Later,  Boucher  (1967)  observed  severe 
storms  to  have  larger  tops  than  non-severe  storms. 

Severe  local  storm  research  becaome  very  popular  in  the  1970s  with 
the  advent  of  geosynchronous  satellites  featuring  ever  improving  spatial 
and  temporal  resolution  in  both  the  visible  and  infrared.  Anvil  bright¬ 
ness  was  examined  by  Arn  (1975)  and  Tuckman  (1982).  Purdom  (1976)  and 
Beckman  (1982)  studied  merging  thunderstorms.  Anvil  growth  rate  has  been 
linked  to  severe  thunderstorms  by  Sikdar  (1970),  Purdom  (1971),  Adler 
and  Fenn  (1976),  Yuen  (1977),  and  Negri  (1977).  Fujita  (1978)  was  able 
to  isolate  an  incident  of  pulsating  anvil  growth,  while  Marshall  and 

Peterson  (1979)  studied  a  storm  whose  anvil  grew  in  area  at  a  rate  of 

2  -1 

700  km  min  yet  did  not  become  severe.  Anderson  (1979)  also  reports 


unusual  anvil  characteristics  while  investigating  anvil  outflow  pat¬ 


terns. 

Adler  and  Fenn  (1979a)  used  short  interval  3-5  minute  infrared 
images  from  geosynchronous  satellites  to  examine  the  rate  of  change  in 
temperature  of  severe  and  non-severe  thunderstorm  tops.  This  led  to  a 
relationship  between  storm  growth  rate  as  compared  to  the  storm’s  up¬ 
draft  velocity  (Adler  and  Fenn,  1979b).  Pryor  (1978)  and  Reynolds 
(1980)  investigated  the  top  temperatures  of  hailstorms. 

The  latest  technique  developed  to  study  thunderstorms  by  satellite 
is  stereographic  observations  (Hasler,  1981}  Fujita,  1982;  Heymsfield 
et  al,  1982).  This  method  features  good  horizontal  and  vertical  resolu¬ 
tion  ( 1  km  x  1/2  km)  but  is  useful  only  in  daylight.  Stereographies 
uses  two  geosynchronous  satellites  focused  on  the  same  cloud  from  dif¬ 
ferent  viewing  positions  to  yield  accurate  measurements  of  height  changes 

History  of  Lidars 

A  lidar  is  a  laser  often  used  for  observing  the  atmosphere.  Lidar 
is  an  acronym  for  light  detection  and  ranging.  The  lidar ’s  principle  of 
operation  closely  parallels  that  of  the  weather  radar  ( radio  detection 
and  ranging) .  Some  researchers  refer  to  the  lidar  as  "laser  radar"  or 
"optical  radar." 

The  lidar  generates  a  short,  powerful  light  pulse  that  propagates 
through  the  atmosphere.  The  light  from  the  laser  pulse  is  backseat tered 
by  atmospheric  gases  and  by  suspended  aerosols  and  is  collected  by  a 
receiving  telescope.  The  telescope  contains  a  narrow  band  interference 


filter  which  removes  background  sky  light  and  directs  a  filtered-returned 
signal  into  a  photomultiplier  tube  where  the  signal  is  enhanced. 

Three  years  after  lasers  became  operational,  Guy  Goyer  and  R.  Watson 
(1963)  are  credited  as  the  first  to  point  out  that  the  laser  is  a  powerful 
tool  for  remote  sensing  of  the  atmosphere.  Shortly  thereafter,  Fiocco  and 
Smullin  (1963)  used  a  ruby  lidar  to  analyze  the  mesosphere.  In  addition, 
Ligda  (1963)  was  the  first  to  use  the  lidar  for  meteorological  purposes  in 
the  troposphere.  Indeed,  the  lidar  has  many  applications  as  a  tool  for 
remote  sensing  of  the  atmosphere  (Collis,  1970;  Derr,  1977).  The  wide 
selection  of  papers  presented  at  the  Eleventh  International  Laser  Radar 
Conference  (NASA,  1982)  emphasizes  the  exceedingly  broad  range  of  potential 
applications  of  lidar  to  the  atmospheric  sciences  from  the  analysis  of 
aerosols  to  wind  measurements. 

NASA  recently  completed  a  survey  of  the  scientific  problems  and  the 
feasibility  of  conducting  lidar  experiments  from  the  space  shuttle.  The 
study  identified  a  number  of  experiments  that  would  be  of  scientific 
interest  and  be  technically  feasible  with  state-of-the-art  technology 
(NASA,  1979). 

According  to  Atlas  and  Korb  (1981)  a  shuttle-borne  lidar  could  con¬ 
duct  experiments  such  as  the  measurement  of  all  atmospheric  state  variables 
(temperature,  pressure,  humidity,  and  wind)  as  a  function  of  altitude 
above  any  surface  location.  Since  this  thesis  examines  a  use  for  a 
shuttle-borne  temperature  sensing  lidar,  this  literature  review  will 


be  restricted  to  laser  methods  of  measuring  temperature. 


Lidar  Sensing  of  Temperature 


Techniques  to  measure  temperature  with  optical  devices  began  with 
Elterman  (1964)  when  he  employed  a  searchlight  to  measure  stratospheric 
density  profiles.  Temperature  profiles  are  computed  from  the  density 
measurements  by  assuming  that  the  atmosphere  is  in  hydrostatic  equilib¬ 
rium.  Twenty-seven  years  later,  Chanin  and  Hauchecorne  (1980)  performed 
essentially  the  same  experiment  using  a  lidar.  Searchlight  procedures 
are  not  accurate  in  the  troposphere  where  all  sensible  weather  occurs, 
but  lidar s  are  quite  useful, 

Raman  lidar  systems  using  rotational  Raman  scattering  from  a  mix  of 
atomspheric  0^  and  gases  cam  recover  temperature  profiles  (Gill  et  al, 
1974;  Cohen  et  al,  1976).  Nevertheless,  Raman  scattering  has  a  low  cross- 
section,  or  probability  of  occurrence.  Hence  this  technique  has  low 
sensitivity  which  makes  it  inappropriate  for  atmospheric  probing  from 
shuttle  altitudes. 

Atmospheric  temperatures  can  also  be  recovered  from  the  three  wave¬ 
length  method.  Such  a  system  transmits  two  wavelengths  centered  on  the 
absorption  lines  of  a  gas  while  the  third  wavelength  is  used  to  calibrate 
the  lidar.  This  method  is  similar  to  the  DIAL  technique  reviewed  below 
(Mason,  1975,  Endemann  and  Byer,  1981).  This  method  is  feasible  from  the 
ground  or  an  airplane  but  not  at  shuttle  altitudes. 

DeLuisi  et  al  (1975)  developed  another  three  wavelength  technique  to 
measure  temperature.  In  addition,  this  method  uses  two  zenith  angles  to 
obtain  separate  density  profiles.  This  technique  is  not  fully  developed, 
and  results  are  very  restrictive,  Russel  and  Norley  (1982)  refined  a  one 


wavelenath 


single-angle  analysis  technique  to  arrive  at  a  density  profile.  This 
system  is  well  suited  for  atmospheric  sounding  from  the  space  shuttle. 

The  final  laser  technique  to  measure  temperature  is  the  DIAL  tech¬ 
nique.  DIAL  is  an  acronymn  for  differential  absorption  lidar .  Korb  and 
Weng  (1979)  theorized  that  one  wavelength  centered  on  a  resonant  absorp¬ 
tion  line  could  be  used  to  make  highly  accurate  temperature  measurements 
and  is  called  the  on-line  wavelength.  A  second  wavelength  is  tuned  to 
an  off-line  frequency  where  there  is  no  molecular  absorption.  The  ab¬ 
sorption  coefficient  at  line  center  is  very  sensitive  to  temperature 
through  the  Boltzmann  distribution.  Oxygen  is  the  ideal  element  to  use 
with  its  wide  absorption  band  near  770  nm.  Furthermore,  oxygen  is  uni¬ 
formly  mixed  throughout  the  atmosphere  (Smith  and  Platt,  1977;  Kolshoven 
et  al,  1980;  Murray  et  al ,  1980;  Lebow  et  al ,  1982) .  This  latter  tech¬ 
nique  yields  very  accurate  temperature  measurements  (less  than  1°C  error) 
and  is  the  technique  employed  in  this  thesis  to  recover  thunderstorm 
cloud  top  temperatures. 


DIAL  Lidar  vs  Infrared  and  Microwave  Radiometers 

The  Global  Atmospheric  Research  Program  (GARP,  1973)  requires  at¬ 
mospheric  temperatures  to  be  +_  1°C  for  2  km  vertical  layers  up  to  the 
tropopause  so  that  numerical  models  forecasting  the  weather  can  be  im¬ 
proved. 

Infrared  and  microwave  radiometers  are  currently  capable  of  1.5  - 
3.0°C  for  a  5-10  km  vertical  layer  (Kolshoven  et  al,  1980).  The  proposed 
new  Advanced  Meteorological  Temperature  Sounder  (AMTS)  will  deliver 
+  1°C  accuracy;  however,  the  vertical  resolution  is  limited  to  4  km. 


These  passive  radiometry  techniques  require  extensive  atmospheric  model¬ 
ing  to  derive  the  temperature  profile  from  the  measured  data.  A  lidar 
in  low  earth  orbit  can  deliver  temperature  errors  less  than  1°C  and 
vertical  resolution  from  10  m  to  2  km  (Atlas  and  Korb,  1981;  Russel, 
1982). 

Scattering  Theory 

There  are  several  scattering  phenomena  that  attenuate  a  laser  beam 
as  it  propagates  through  the  atmosphere.  These  phenomena  are  Rayleigh 
scattering,  Mie  scattering,  resonant  absorption,  and  Raman  scattering 
(Ely,  1972).  This  thesis  will  consider  the  first  three  scattering 
phenomena  which  are  tabulated  at  different  wavelengths  and  altitudes  by 
McClatchey  et  al  (1971). 


Rayleigh  Scattering 

Rayleigh  theory  considers  the  scattering  of  incident  radiation  off 
particles  whose  diameters  are  very  much  smaller  than  the  incident  wave¬ 
length.  The  irradiance  of  light  can  be  measured  at  a  distance  r  from  a 
scattering  atom  or  molecule  and  expressed  as 


where  e  =  irradiance  of  incident  radiation  (watt  m"2) 
o 

=  Rayleigh  cross  section  (m^). 


(1) 


Kerker  (1969)  explains  that  Rayleigh  theory  represents  the  irradiance 


(2) 


e  =  e  2ir  (n-1)  (1  +  cos  0} 
o _ 

2  »,2  i4 

r  N  A 

where  X  =  wavelength  of  incident  light 
N  =  particle  number  density 
n  =  apparent  index  of  refraction  of  the  gas 
0  =  scattering  angle  in  the  forward  direction. 


Thus,  the  Rayleigh  cross  section  per  solid  angle  is 


dar  =  2tt2  (n-1)2  (1  +  cos20)  (3) 

“dO  N2  X4 


or 


2ir  sin  0  d  0. 


(4) 


Mie  Scattering 

Mie  theory  considers  the  scattering  of  incident  radiation  off 
particles  about  the  same  size  as  the  incident  wavelength.  Expressions 
for  Mie  theory  are  complicated,  since  the  scattered  irradiance  is  a 
function  of  scattering  angle,  ratio  of  particle  diameter  to  wavelength, 
and  the  particle's  complex  index  of  refraction  (Ely,  1972).  The  ir¬ 
radiance  of  Mie  scattered  light  at  a  distance  r  from  a  scatterer  can  be 
written  as 


e 


*2  ii> 


2  2 
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(5) 


where  e  =  irradiance  of  incident  radiation 
o 

A  =  wavelength  of  incident  radiation 

i^  and  i^  =  irradiance  functions  perpendicular  and  parallel  to  the 
scattering  plane. 


Equation  1  can  be  rewritten  for  Mie  theory  as 


e  =  £o  gm  (6) 

2 

r 


where  0  =  Mie  theory  cross  section.  Thus,  0  can  be  written  as 

m  m 


a 

m 
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(7) 


McClatchey  et  al  (1971)  points  out  that  the  Mie  cross  section  0 
—  —  m 

consists  of  the  sum  of  an  absorption  cross  section  and  a  scattering 

cross  section  0  such  that 
s 


0=0+0.  (8) 

mas 

McClatchey  also  computes  extensive  tables  using  Mie  theory  when  computing 
aerosol  extinction  at  most  laser  wavelengths  using  different  model 
atmospheres  and  altitudes. 


Resonant  Absorption 

Resonant  absorption  is  generally  referred  to  as  the  gaseous  equiva¬ 


lent  to  Mie  theory.  In  Mie  theory,  the  imaginary  component  of  the  index 


of  refraction  accounts  for  0  ;  however,  in  resonant  theory  absorption  is 


characterized  by  an  absorption  coefficient  a  .  Here  the  first  subscript 

gv 

g  refers  to  an  absorbing  gas  and  the  second  subscript  refers  to  the  par¬ 
ticular  frequency  at  which  absorption  is  taking  place. 

Every  molecule  in  the  atmosphere  has  certain  allowed  transitions  bet 
ween  its  electronic,  vibrational,  and  rotational  energy  states.  Accord¬ 
ingly,  the  absorption  coefficient  can  be  very  large  at  certain  frequency 
bandwidths . 

A  laser  beam  propagating  through  the  atmosphere  will  have  photons 
absorbed  by  atmospheric  molecules  and  removed  from  the  beam  if  the 
photon’s  energy  equals  that  of  a  molecular  energy  transition.  In  equa¬ 
tion  form 

E  =  hv  (9) 

where  E  =  energy  needed  for  molecular  transition 
hv  =  energy  of  a  photon  from  laser  beam. 

The  many  possible  transitions  of  particular  molecule  are  referred  to 
as  lines  or  resonances.  These  lines  are  not  infinitely  sharp  in  a  freq¬ 
uency  spectra,  but  are  collision  broadened  for  altitudes  up  to  the  trop- 
opause. 

There  is  an  extensive  development  of  resonant  absorption  a ,  in 

Chapter  IV  of  this  thesis. 


III.  Lidar  Equation 


A  lidar  mounted  in  the  cargo  bay  of  the  space  shuttle  can  transmit 
power  earthward  and  incident  on  a  volume  located  at  the  top  of  a  thunder¬ 
storm.  The  power  arriving  at  the  cloud  top  can  be  expressed  by  Bouger's 
law  as 


P  (z)  =  P  exp 
T  o 


dz 


where  Z  =  height  of  shuttle  above  surface  ()un) 

=  height  of  thunderstorm  top  (km) 
p^  =  transmitted  power  (watts) 

P^  =  power  at  thunderstorm  top  (watts) 
a(z)  =  total  atmospheric  extinction  coefficient  (km  ^) 
ja(z)dz  =  atmospheric  attenuation. 


(10) 


Some  of  this  laser  power  reaching  the  cloud  top  will  be  backscattered 

to  space.  This  backscattered  power  is  the  product  of  the  backscatter 

coefficient  for  the  volume  being  sensed  and  the  depth  of  this  volume. 

The  total  power  per  solid  angle  backscattered  to  space,  which  is  the 

intensity  I  ,  can  be  expressed  as 
B 


I  =  p_(z)  B(z)d 


(11) 
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where  IB  =  intensity  backscattered  to  space  (watts  sr 
P^tz)  =  power  at  thunderstorm  top  (watts) 

B(z)  =  volume  backscatter  coefficient  (km  ^sr 
d  =  vertical  resolution  (depth  of  volume)  (km) 
C  =  speed  of  light  (km/sec) 
tp  =  pulse  width  (sec) . 


The  amount  of  laser  backscattered  power  that  actually  enters  the 
lidar's  receiving  telescope  can  be  expressed  as 

z 

Pr (z )  =  ^r  In  exp  -  f  a(z)dz  (13) 

R2  Jz 

T 

where  Pl(z)  =  power  received  by  lidar  (watts) 

Ar  =  area  of  receiving  telescope  (km) 

R  =  distance  from  cloud  top  to  lidar  (km) . 


Ar 

Note  that  — ~  is  the  solid  angle  subcended  by  the  receiver  at  range  R, 
since  lidar  systems  are  usually  configured  with  the  divergence  angle  of 


the  transmitted  beam  less  than  the  receiver's  field  of  view. 


All  of  these  equations  can  be  combined  to  yield  the  lidar  equation 


for  a  shuttle-borne  lidar 


PL(z)  =  PQ  e  Ar  Ctp  B (z)  exp 
2  R2 


(14) 


where  e  is  the  lidar's  system  efficiency,  the  product  of  the  optical  sys¬ 
tem  efficiency  and  the  receiving  system  efficiency.  The  volume  backscatter 
coefficient  B(z)  is  defined  as  the  fractional  amount  of  incident  power 


scattered  per  steradian  in  the  backward  direction  per  unit  path  length. 


The  volume  backscatter  coefficient  is  the  sum  of  scattering  contributions 


from  all  atmospheric  constituents  (gases,  aerosols,  and  cloud  ice  crystals) 
and  includes  both  inelastic  and  elastic  processes.  It  can  be  expressed  as 

n  n 

B(z) =  Z  Bi<z)  =  I  N±  dai(z)  (15) 

i=l  i=l  — 3ft — 

where  i  =  summation  over  all  scattering  processes 

Nj_  =  number  density  involved  in  ifch  process  (km-2) 
da(z)/dft  =  backscattering  cross  section  of  the  ifch  process  (km2sr-^) 

B(z)p  =  backscatter  coefficient  (km"^sr"l) . 

So  in  equation  14,  B(z)  can  be  replaced  with  Bg  +  BA  +  Bc  where  the  sub¬ 
scripts  apply  to  gases,  aerosols,  and  cloud  particles  respectively. 

The  total  extinction  a(z)  is  the  sum  of  the  extinction  coefficients 
of  all  processes  in  the  atmosphere  hence 

a(z)  =  cig(z)  +  otA(z)  +  ac(z)  (16) 

where  ag  =  extinction  due  to  atmospheric  gases  (resonant  absorption)  (km_^) 
aA  =  extinction  due  to  aerosols  (km-l) 
ac  =  extinction  due  to  cloud  particles  (km-^) . 

Furthermore,  each  of  these  coefficients  is  composed  of  a  contribution 
from  scattering  and  absorption.  The  equation  for  total  extinction  can  be 
rewritten  as 

a  =  Sg  +  Kg  +  SA  +  Ka  +  Sc  +  Kc  (17) 

where  =  scattering  due  to  the  ifch  constituent  (km“^) 

Ki  =  absorption  due  to  the  constituent  (kr.»“-) . 


and  i  is  either  for  gases,  aerosols,  or  cloud  particles  (g,  A,  c) . 

Combining  the  above  results  with  equation  14,  the  lidar  equation  for 
an  instrument  mounted  on  the  shuttle  can  be  expressed  as 


PL(z)  -  P0  e  Ar  C  tr 


Bg(z)  +  BA(z)  +  Bc(z) 


exp 


fZL 

-2  [ctg(z)+aA(z) 
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+  ac (z)  ]  dz 


(18) 


Multiple  Scattering 

The  lidar  equation  expressed  by  equation  18  considers  the  theory  of 
single  scattering  as  the  beam  scatters  off  gases,  aerosols,  and  cloud 
particles.  Actually  equation  18  is  only  an  approximation  since  multiple 
scattering  prevails  in  reality.  Multiple  scattering  greatly  increases 
the  backscattered  signal  as  the  laser  beam  enters  a  dense  scattering  situa¬ 
tion  such  as  a  cloud  top. 

Figure  6  illustrates  the  multiple  scattering  phenomena.  From  shuttle 
altitude  and  down  to  a  cloud  top,  weak  single  scattering  prevails  as  the 
laser  illuminated  volume  is  entirely  restricted  to  the  volume  enclosed 
within  the  divergence  angle  9  of  the  laser  beam.  In  other  words,  single 
scattering  occurs  when  G  is  less  than  the  receiver's  field  of  view  <p. 
Multiple  scattering  occurs  as  the  beam  enters  the  cloud  top  and  spreads 
as  shown  in  Figure  6  by  the  shaded  region.  Photons  in  this  shaded  region 
require  three  or  more  scatterings  before  being  directed  back  to  the  re¬ 
ceiver  since  this  region  is  greater  than  <j>  (Hobbs  and  Deepak,  1981)  . 

Multiple  scattering  is  neglected  in  this  thesis  and  equation  18  will 


not  be  modified. 


IV.  Dial  Theor 


The  laser  technique  used  to  measure  temperature  in  this  thesis  is 
the  differential  absorption  lidar  (DIAL)  technique  (see  literature  review) 
A  lidar  on-board  the  space  shuttle  using  the  DIAL  technique  will  transmit 
two  wavelengths  along  the  same  path  down  to  the  top  of  a  thunderstorm. 

The  first  wavelength,  an  on-line  wavelength,  is  selected  to  match  the 
center  of  a  resonant  absorption  oxygen  line  in  the  oxygen  "A"  band  near 
770  nm.  The  second  wavelength,  an  off-line  wavelength,  is  used  as  a 
reference  and  is  chosen  at  a  nearly  spectral  position  where  there  is 
almost  no  absorption  from  the  oxygen  "A"  band.  This  DIAL  technique  is 
designed  so  that  scattering  and  absorption  properties  of  the  atmosphere 
are  identical  except  for  the  resonant  absorption  at  the  on-line  wave¬ 
length. 

Remote  Sensing  of  Temperature 

Kalshoven  et  al  (1980)  outlined  a  procedure  to  determine  atmospheric 
temperature  measurements  by  using  the  DIAL  technique.  This  technique 
will  be  modified  by  the  author  so  that  thunderstorm  cloud  top  temperatures 
can  be  recovered  by  operating  a  DIAL  lidar  from  the  space  shuttle . 

The  differential  absorption  coefficient  (extinction  due  to  gases) 

Gig  can  be  found  by  the  ratio  of  power  on-line  to  power  off-line  returned 
to  the  lidar' s  receiving  telescope.  This  can  be  expressed  in  equation 
form  developed  by  Kalshoven  et_  al  (1980)  and  modified  by  the  author  as 


where  z  =  distance  between  lidar  and  cloud  top  (km) 

6.9  =  scale  height  for  an  isothermal  atmosphere  (km) 


cig0  =  differential  absorption  coefficient  at  line  center  or  extinction 
due  to  gases  (km-^) 


P  =  received  on-line  power  (watts) 


P  =  received  reference  power  (watts). 


The  differential  absorption  coefficient  hereafter  called  the  ab¬ 
sorption  coefficient  or  extinction  coefficient  due  to  gases  can  be  found 
analytically  with  an  expression  derived  by  Elasser  (1960)  as 
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where  aQ( v)  =  absorption  coefficient  at  wave  number  v 


S  =  line  strength 


a  =  collisional  broadening  half-width  at  half  maximum  at 


standard  pressure  and  temperature  or  273.16°k,  1013.25  mb 


Ps  =  pressure  at  STP 


P  =  pressure  at  specified  level 


=  temperature  at  STP 


T  =  temperature  at  specified  level 
v— v0  =  frequency  bandwidth  of  absorption  line 
vo  =  wave  number  of  center  of  absorption  line 


STP  refers  to  conditions  at  standard  temperature  and  pressure. 


A  more  conventional  form  of  this  equation  for  use  in  the  atmosphere 


where  pressure  and  temperature  change  with  height  is 


C(go  =  N(P,T)  S  (T)  F  (\T-\J0) 


(21) 


where  N(P,T)  =  atmospheric  number  density  for  oxygen  (molec  m”3) 

S  (T)  =  line  strength  (rnolec-^  m2  cm-3-) 

P(^v0)  =  line  shape  function  (cm). 

The  number  density  can  be  found  from 

N(P.T)  =  q  Wns  f-  jjjS-  (22) 

where  Ps  =  pressure  at  STP  (1013.25  mb) 

P  =  pressure  at  specified  height  (mb) 

Ts  =  temperature  at  STP  (273.16ok) 

T  =  temperature  at  specified  height  (°K) 
q  =  dry  air  mixing  ratio  for  O2  (0.209) 
ns  =  dry  air  number  density  at  STP  (2.69  x  lO2^  molec/m3) 

W  =  correction  factor  for  atmospheric  water  vapor. 

Kalshoven  et^  a]L  (1980)  was  able  to  derive  an  empirical  expression  for  W 
given  as 

W  =  1  -  2.23  x  109  ~  exp  (  -  52®^  (23) 

where  RH  =  relative  humidity  (.80  for  80%). 

W  was  found  to  be  nearly  one  and  contributes  <  2%  to  the  computation  for 
number  density. 

Line  strength  S  for  oxygen  can  be  found  from  Penner  (1959)  as 


cig(v)  dv 


(24) 
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where  ve=  wave  number  of  resonant  line  center 
M  =  mass  of  the  molecule 
Q(T)  =  the  partition  function 
T  =  temperature 

E"  =  energy  in  wave  numbers  of  the  lower  state 
J"  =  total  angular  momentum  quantum  number 
f(J")  =  statistical  weight 

Rev  *  combined  electronic-vibrational  transition  moment 
c  =  speed  of  light 
h  =  Planck's  constant 
k  =  Boltzmann's  constant 
*  =  multiplication  operator. 


Burch  and  Gryvnak  (1969)  were  able  to  simplify  this  expression  and  derive 
an  expression  for  line  strength  of  the  oxygen  "A"  band  at  a  particular 
line  as 


Sj» 


Vj" 

Q  (T) 


exp  (-he  E"/k  T) . 


(26) 


The  rotational  partition  function  Q(T)  is  given  by 


Q(T)  =  £  Fj"  exp  (-h  c  E"/k  T) 


(27) 


where  Fj«  is  a  weighting  function. 


Finally,  McClatchey  et  al_  (1971)  expressed  line  strength  in  a  very  opera¬ 
tional  format  as 


Ts 


S(T)  =  S  (Ts)  —  exp 


1.439  E"  (  ~~  -  i) 
Ts  T 


(28) 


where  S(Ts)  =  O2  line  strength  at  STP  (molec-^  m2  cm-*) 

E"  =  energy  in  lower  level  of  the  oxygen  molecule  (cm~l) . 


Burch  and  Grynak  (1969)  have  compiled  extensive  tables  of  S(Ts),  E”, 
J”,  etc.  for  lines  throughout  the  oxygen  ”A”  band. 


The  line  shape  function  F(v-v0)  is  a  combination  of  Lorentz  and  Doppler 
broadening  and  is  typically  expressed  by  the  Voigt  integral  as 
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where  x  =  (In  2 )^2 
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a  =  — ^  (In  2)  l/2 
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and  aD  =  Doppler  broadened  halfwidth ( cm  ) 


_1 


=  Lorentz  broadened  halfwidth(cm  ") 


ressions  for  ap  and  aj,  are 
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where  m  =  molecular  mass  for  O2  (5.31  x  10“26  kg) 
a®  =  collision  halfwidth  for  O2  at  STP. 


According  to  Penner  (1959)  and  Kalshoven  et  al  (1980),  the  line  shape 


function  F(v-v0)  reduces  to 


F  (0)  =  F '  exp 


[K 


exp  (-2^) dz 


where  F(0)  =  line  shape  at  line  center 
z  =  height 


a  =  correction  factor. 


The  value  of  the  correction  factor  "a"  varies  with  the  molecule  of 


interest.  For  oxygen 


c  ntS  p 

a  =  2.14  x  105  L 


where  is  expressed  in  cm“l 
P  is  expressed  in  mb 
T  is  expressed  in  °K 
vD  is  expressed  in  cm“l. 


Solving  the  Voigt  integral  can  be  avoided  by  using  a  good  analytical 


approximation  such  as 


At  heights  of  typical  thunderstorms  this  equation  is  93%  accurate 
(Penner,  1959). 

Recalling  Equation  21,  and  combining  all  of  the  above  equations, 
the  absorption  coefficient  is  calcinated  by 


og0  =  N(P,T)  S(T)  F( v  -  vQ) 


n"3/2,P  ,  ,Ts 


Is, 


=  q  no  W  T  '"‘(— )  (-^)  S(Ts)  l-f)  exp 


1.439  E" 

Ts  T 


1 

3ag 


n  _ 

1  10a 


1  1  P  Ts  *5 

where  4s  =  a -  <p->  <^r) 

L  L  S 


(38) 

(39) 


Therefore, 


-3/2 

ago  =  q  ns  w  T  SIT  )  exp 


1  439  E"  (— —  -  — )  -i—  (  1  _  exP  ^ 

1.4J9  E  'ts  T  ^  (1  10a 


7x5 

L 


—3/2  F" 

Aq  T  '  W  exp  (-1.439 


(1  - 

1  10a  J 


(40) 


where  Aq  =  T_3/2  -S(4§}  exp  (1.439  ~)  .  °K3/2  m"1  (41) 

L  S 


This  equation  can  be  solved  for  temperature  and  arrive  at  Kalshoven’s 
et  al  (1980)  iterative  expression  for  temperature.  Equation  40  can  be  re¬ 
written  as  an  iterative  expression 


“go 


=  A,  T,-3/2  W  exp  (-1.439  E"/ti+1>  (1  -  ~^(-1>) 


o  1 


10a 


(42) 


Now  solving  for  T^+^  so  that  an  equation  yielding  a  laser  derived  temp¬ 


erature  can  be  founds 


Where  is  a  seed  temperature  value  (a  standard  atmosphere  value)  and  is 
not  critical  to  the  results  of  Ti+l.  Kalshoven  et  al^  (1981)  tested  this 
DIAL  technique  at  sea  level  using  the  .7684  um  O2  line  and  found  a 
temperature  accuracy  of  0.6°K. 

To  investigate  the  rate  of  change  of  absorption  with  respect  to 
temperature,  the  first  derivative  of  equation  40  can  be  taken  with  re¬ 
spect  to  temperature.  Equation  40  is  rewritten  for  convenience 

“g0  =  Aot3/2  w  exp  (“  1-439  E"/T) 

It  is  appropriate  to  omit  the  weak  pressure  components  to  0igo  such  that 


1  exp  (1) 
10a 


m-3/2 


,-1.439  E"  , 


Now  taking  the  derivative 


dagQ  =  clT  ( 

-  T 

ago 


1.439E" 

T 


(47) 


This  equation  says  that  the  shuttle  lidar  sensitivity  to  temperature 
changes  improves  by  choosing  higher  energy  lines  within  the  oxygen  "A" 
band.  However,  if  E"  is  increased  by  choosing  another  line,  the  line 
strength  becomes  far  too  weak,  and  the  absorption  coefficient  becomes 
too  small.  A  trade-off  must  be  made. 

The  on-line  wavelength  chosen  in  this  thesis  is  .7696  um.  This 
choice  lies  within  that  portion  of  the  oxygen  "A"  band  offering  optimum 
temperature  sensitivity  and  absorption  (Kalshoven  et  al ,  1980) . 

The  off-line  or  reference  wavelength  is  chosen  to  be  .7614  um  since 
oxygen  absorption  at  this  wavelength  is  negligible  compared  with  that  at 
the  on-line  wavelength  (Smith  and  Platt,  1977) . 


Relative  Humidity 

Relative  humidity  RH  or  atmospheric  water  vapor  will  slightly  affect 
the  calculation  for  the  number  densigy  of  O2.  RH  is  defined  as  the 
ratio  of  the  actual  mixing  ratio  of  a  parcel  of  air  at  a  specific  pres¬ 
sure  and  temperature  to  the  mixing  ratio  that  saturated  air  would  have 
at  the  same  pressure  and  temperature  or 

RH  ■  -  (100)  ’  (48) 

ws 

where  w,  ws  are  the  mixing  and  saturated  mixing  ratio  respectively.  A 
mixing  ratio  is  merely  the  mass  of  water  vapor  contained  by  a  mass  of 
dry  air  (Hess,  1956) . 

The  calculation  of  the  O2  number  density  N(P,T)  (see  Equation  22) 
requires  a  correction  for  atmospheric  water  vapor  expressed  in  Equation 
23.  The  values  of  RH  required  in  this  equation  will  vary  according  to 
height  (AWSM  105-124,  1969)  as 


Surface  -  3 

km 

50% 

4-6 

km 

30% 

O 

H 

1 

r* 

km 

10% 

11  -  20 

km 

5% 

Differential  Absorption  Coefficient  Calculation 

The  differential  absorption  coefficient  for  molecular  oxygen  at 
.7696  pm  can  now  be  calculated  using  the  previously  developed  DIAL 
theory  in  a  sequential  manner  consisting  of  seven  steps. 

This  example  will  calculate  the  absorption  coefficient  OgQ  at  10  km 


altitude  and  will  demonstrate  the  appropriate  units  to  use. 


The  following  spectral  information  is  provided: 


X  =  .7696  ym 
vQ  =  12,988.734  cm-1  * 

E"  =  1421.436  cm-1  * 

S(Ts)  =  3.377  x  10-30  molec"3-  m2  cm“l  * 

=  .038  cm-*-  * 

Ps  =  1013.25  mb 

P  =  281  mb  (McClatchey 1 s  et  al  (1971)  mid-latitude  summer  atmosphere) 
Ts  =  273 . 16°K 

T  =  235°K  (McClatchey ' s  et  al^  (1971)  mid-latitude  summer  atmosphere) 

where  *  indicates  values  tabulated  by  Burch  and  Gryvnak  (1969) . 

The  resonant  absorption  coefficient  oig0  or  extinction  due  to  oxygen 
"A”  band  absorption  ag  (km-1)  (subscription  g  refers  to  02  gas)  are 
names  for  the  same  quantities.  Recall  that  agQ  can  be  calculated  from 
Equation  21  as 

a9o  =  N(P,T)  s  (T)  F(V-V0) 

where  the  following  seven  steps  are  required  to  solve  the  expression. 

A.  A  correction  for  atmospheric  water  vapor  W  must  be  found  from  Equation 
23. 

=  .9999 


W  =  1  - 


2.23  x  10- 


.10 


281  mb 


exp  (- 


5385 

285°k 


B.  The  number  density  of  O2  at  10  km  can  now  be  calculated  from  Equation 


N(P,T)  =  (.209)  (.9999)  (2.69xl025  molec  m-3)  jg) 


24  -3 

=  1.8126x10  molec  m 


C.  Line  strength  S(T)  using  Equation  28 

S(T)  =  (3o377xl0~30  molec”*  m2  cm”*)  ( ^ )  exp 

,  1  _  _ 

1 273 . 16°K  235°K' 

=  1.5082x10  3^  molec  *  m2  cm  * 


1.439(1421.436) 


D.  The  Lorentz  broadened  half width  a  for  the  resonant  line  using 
Equation  34. 


o  =  (0.038  cm 

Xj 


-1 


)  ( 


281  mb 
1013.25  mb 


,  ,  273 , 16°K,  1/2  .  -1 

)  (~235ck  =  °-0114  cm 


E.  The  weak  pressure  correction  factor  "a"  using  Equation  36, 

„  .  „  ,„5  ,  0.038  cm-1  ,  .281  mb,  „ 

a  -  2.14x10  ( 12j988o734  Cm-1)  1 235°K  }  =  0,7486 

F.  Line  shape  at  line  center  F(0)  using  Equation  37. 


F(v 


-vQ)  -  f(o)  =  (3)  (oToirrHiFTy 


1  - 


exp  (1) 


(10)  (0.7486) 


=  17.4012  cm 


G.  The  absorption  coefficient  dg0  can  now  be  calculated  using  Equation  21. 


24  -3  -30  -1  2  -1 

ag0  =  (1.8126x10  molec  m  )  (1.5082x10  molec  m  cm  ) 

(17.4012  cm)  ( 103  m  km  *)  =  ,0513  km”* 

This  sequential  scheme  of  seven  steps  can  readily  be  programmed  using 


a  suitable  computer  language.  The  absorption  coefficient  was  computed  in 
this  manner  from  the  surface  to  17  km  using  McClatchey's  et  al  (1971) 


mid-latitude  summer  atmosphere .  These  values  are  displayed  in  Table  I 
which  compare  to  Kalshoven's  et:  al_  <1980)  work  when  extended  to  17  km  and 
using  an  on-line  wavelength  of  .7684  ym.  Figure  7  displays  agQ  as  a 
function  of  height. 

Predicting  with  Polynomial  Regression 

It  may  be  of  interest  to  calculate  the  02  extinction  coefficient 
using  regression  analysis  to  statistically  model  the  relationship  between 
the  independent  variable  (height)  and  a  dependent  variable  (ctgQ)  . 

The  exact  form  of  the  true  functional  relationship  between  the  02 
extinction  coefficient  and  height  is  assumed  to  be  unknown  so  it  is 
appropriate  to  use  a  polynomial  model  as  the  approximating  function. 
Regression  analysis  of  this  nature  assumes  that  height  can  be  measured 
with  negligible  error  {Hines  and  Montgomery,  1980) . 

To  determine  the  relationship  between  the  two  variables  a  cubic 
equation  will  be  fitted  using  the  least  squares  method.  This  cubic 
equation  can  be  written  as 

ago  =  a  +  bz  +  cz2  +  dz2  (49) 


V.’  *. 


where  z  =  height  (km) 

ag0  =  O2  absorption  coefficient  or  extinction  due  to  gas  at  line 
center  (km-1). 


The  regression  coefficients  (a,b,c,d)  can  be  found  by  solving  the 
following  system  of  equations.  These  equations  can  be  solved  numerically 
using  Gaussian  elimination  with  partial  pivoting  (Ketter  and  Prawel,  1969): 


i 
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Table  I 


Oxygen  Absorption  Coefficient  km  )  from  Surface  to  17  km 


altitude  (km) 


a  ( km  ^ ) 

_ 


.2848 


.2604 


.2326 


.2023 


.1746 


.1488 


.1254 


.1042 


.0832 


.0669 


.0513 


.0393 


.0284 


.0202 


.0164 


.0118 


.0066 


.0005 
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V.  Atmospheric  Consideraions 

Index  of  Refraction 

Edlen  (1953)  computed  the  index  of  refraction  at  1013  mb  and  288°k 
for  any  particular  wavelength  as 


(n-l)10-8 


6432.8  + 


2949810  25540 

146- ( a-^ )  +  41- ( A- ^ ) 


(53) 


where  n  =  index  of  refraction 
X  =  wavelength  (ym) . 


Likewise,  Barrell  and  Sears  (1939)  found  the  index  of  refraction 
of  moist  normal  air  in  the  range  10  to  30°C  at  960  to  1066  mb  for  any 
particular  wavelength  as 


(n-l)lO6  =  0.378125  + 


0.0021414_+  0.00001793 


*fp  1  +  (1-049  -  0. 
[  1  +  0. 

-[ 


F  T  W 

0157  T)P  10-6 


003661  T 


„  0.000680 

0.0624  -  - - 


1  +  0.003661  T 


(54) 


where  T  =  air  temperature  (°c) 

P  =  barometric  pressure  (mm  Hg) 
f  =  vapor  pressure  (mm  Hg) . 


Elterman  (1968)  computed  the  index  of  refraction  n  at  various  wave¬ 
lengths  using  Equation  53.  He  safely  neglects  the  effects  of  water  vapor 
by  working  in  and  around  the  visible  spectrum.  Elterman's  (1968)  value 
for  n  at  .8  ym  will  be  used  in  this  thesis  for  calculations  at  .7696  and 


.7614  ym.  This  value  is 


n{.8  pm)  =  1.00027503. 


Rayleigh  Cross  Section 

Rayleigh  theory  (see  literature  review)  applies  when  laser  radiation 
scatters  off  atmospheric  gases.  The  probability  of  scattering  off  an 
individual  gaseous  constituent  is  the  cross  section.  The  Rayleigh 
cross  section  ar  can  be  expressed  by 


crr(A)  = 


8tt3  (n-1)  2  6  +  3<S 

3X4  Ns2  6-76 


where  cr  =  Rayleigh  scattering  cross  section  (cm2) 

Ns  =  molecular  number  density  at  STP  (cm~3) 

X  =  wavelength  (cm) 

6  =  depolarization  factor  (.035  from  Gucker  and  Basu,  1953) 


Elteman  (1968)  reports  that  the  expression  6  +  36/6  -  76  accounts  for 
the  depolarization  caused  by  the  anisotropy  of  atmospheric  molecules. 
He  further  calculates  or  at  .8  pm  which  will  be  used  at  the  on  and  off 
line  wavelengths.  This  value  is 


Or  (.7696  pm,  .7614  pm)  =  9.990xl0~28  cm2. 


Atmospheric  Extinction 

The  total  extinction  a  expressed  in  Equation  16  consists  of  three 
components.  The  atmospheric  extinction  component  ctg  is  due  to  scattering 
and  absorption  by  gas  molecules  and  is  discussed  in  this  section.  This 
extinction  due  to  atmospheric  gases  oig  can  be  calculated  at  different 
heights  according  to 


CJg  =  or ( A)  N  (z)  (10^  cm  km“l)  (56) 

where  ag  =  extinction  due  to  gas  (km-^) 

N  =  molecular  number  density  which  is  a  function  of  height. 

The  reader  should  note  that  the  above  computation  for  ag  is  correct 
at  the  nonabsorbing  wavelength  (.7614  ym)  .  At  the  absorbing  wavelength 
(.7696  yin),  the  O2  absorption  coefficient  ago(z)  (see  the  remote  sensing 
of  temperature  section)  must  be  added  to  ag  to  obtain  the  extinction  at 
an  absorbing  wavelength.  Values  for  the  Rayleigh  molecular  number  density 
N  can  be  obtained  from  the  U.S.  standard  atmosphere  or  McClatchey  et  al 
(1971)  . 

Table  II  lists  the  molecular  number  density,  atmospheric  extinction 
coefficient  at  the  off-line  wavelength  (.7614  ym) ,  and  the  extinction 
coefficient  at  the  on-line  wavelength  (.7696  ym)  where  O2  absorption 
values  from  Table  I  have  been  added  to  the  computation  for  ag.  The 
heights  ranging  from  0-17  km  were  arbitrarily  chosen;  however,  this  thesis 
is  only  concerned  with  values  from  10  -  17  km  which  are  average  top 
heights  of  severe  and  non-severe  thunderstorms .  It  is  possible  for  thunder 
storms  to  exceed  17  km.  If  ago  is  needed  for  such  rare  storms  than  a  new 
O2  line  should  be  chosen  so  that  agD  is  larger  and  therefore  measureable 
at  extreme  altitudes.  An  ideal  resonant  absorption  line  in  this  case 
would  be  the  .7677  ym  (13,021.283  cm“l)  oxygen  line.  Once  again  a  trade¬ 
off  is  necessary.  If  .7677  ym  is  used  to  measure  super  high  thunderstorms, 
then  system  signal-to-noise  ratio  (discussed  in  Chapter  VII)  becomes  use¬ 
less  for  smaller  storms  because  of  extreme  absorption  at  lower  altitudes 


from  O2. 


Table  II 


Atmospheric  Number  Density  N,  Off-line  (.7614  pm),  and  On-line  (.7696  pm) 
Extinction  due  to  Gas  ag  all  as  a  Function  of  Height  z 


Height  z 

N 

cig  (.7614  pm) 

(.7696  pm) 

(km) 

(cm- 3) 

km-1 

km-3- 

0 

2 . 547xl019 

2.544xl0-3 

.2873 

1 

2.311 

2.309 

.2627 

2 

2.093 

2.091 

.2347 

3 

1.891 

1.889 

.2042 

4 

1.704 

1.702 

.1763 

5 

1.531 

1.530 

.1503 

6 

1.373 

1.371 

.1268 

7 

1.227 

1.225 

.1054 

8 

1.093 

1.092 

.0843 

9 

9.712xl018 

9.702xl0-4 

.0679 

10 

8.598 

8.589 

.0522 

11 

7.585 

7.577 

.0401 

12 

6.486 

6.478 

.0290 

13 

5.543 

5.538 

.0203 

14 

4.738 

4.733 

.0169 

15 

4.049 

4.045 

.0122 

16 

3.461 

3.458 

.0069 

17 

2.959 

2.956 

.0007 
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The  volume  backscatter  coefficient  Bg  (km-^  sr-3-)  for  atmospheric 
gases  in  the  lower  100  km  of  atmosphere  can  be  found  according  to 
Hinkley  (1976)  from 


Bg  =  Br  =  N  5.45 


A  (ym) 
.55 


4.09 


xlO-23 


(57) 


where  Br  =  Rayleigh  b  ackscatter  coefficient  (km~^sr~l) 
N  =  atmospheric  number  density  (cm~3) 


A  =  optical  wavelength  (pm)  . 


Values  for  N  can  be  taken  from  Table  II  or  calculated  by  using  a 
suitable  model  atmosphere  such  as  McClatchey's  et  al  (1971)  mid-latitude 
summer  atmosphere  reproduced  in  part  in  Table  III.  N  can  be  calculated 
by 


N 


pNo 

A 


(58) 


where  p  =  atmospheric  density  which  is  a  function  of  pressure  and 
temperature  (g  cm-3) 

N0  =  Avogadro's  number  (6.023xl023) 

A  =  gram  atomic  weight  of  dry  air  (28.699  g) . 


The  exponent  on  Equation  57  is  4.09  instead  of  4.0  as  prescribed 
by  Rayleigh's  theory  (see  Equation  2)  to  account  for  the  slight  wave¬ 
length  dependence  of  the  refractive  index  of  air. 

Equation  57  is  appropriate  in  all  cases  except  at  those  wavelengths 


where  molecules  have  absorption  coefficients  much  greater  than  the  oxygen 
"A"  band  (Hinkley,  1976) . 


Using  Equations  57  and  58,  Bg  is  found  using  McClatchey's  et  al 
(1971)  mid-latitude  summer  atmosphere.  Such  values  are  presented  at  both 
the  on-line  wavelengths  and  the  off-line  wavelengths  in  Table  IV. 

Attenuation  by  Gases 

In  this  thesis,  a  highly  coherent  laser  beam  is  transmitted  from 
an  orbiting  space  shuttle  down  through  the  atmosphere.  Unfortunately, 
the  beam  will  be  attenuated  every  kilometer  along  its  path  as  photons 
are  scattered  and  absorbed  or  otherwise  removed  from  the  beam.  The 
following  sections  will  develop  relationships  in  an  attempt  to  model  laser 
beam  attenuation  with  height  caused  by  atmospheric  gases,  aerosols, 
and  cloud  particles. 

The  eventual  development  of  an  equation  expressing  how  gases  at¬ 
tenuate  a  laser  beam  calls  for  an  initial  assumption.  It  will  be  assumed 
that  the  atmosphere  can  be  modeled  as  isothermal  from  10  km  to  the  "top" 
of  the  atmosphere  (about  80  km) . 

Using  McClatchey's  et  al  (1971)  mid-latitude  summer  atmosphere, 
temperature  at  10  km, T^q  =  235°k  while  temperature  at  80  km,  Tgg  =  2  16°k. 
So  the  variation  of  temperature  with  height  is 

.919  <  <  1 

T10 

Thus,  an  isothermal  assumption  is  somewhat  reasonable.  According  to  Hess 
(1959) ,  pressure  decreases  exponentially  with  height  in  an  isothermal 
atmosphere  by 

P  =  Po  exp  (  -  £) 


(59) 


Table  IV 


Volume  Backscatter  Coefficient  Bg  (km~l  sr' 
for  Mid-Latitude  Summer  Atmosphere 


Height  On-Line 


km  (km) 

.7696  urn 

0 

3.416xl0"4 

1 

3.097xl0"4 

2 

2.798xl0"4 

3 

2.537xl0"4 

4 

2.294xl0"4 

5 

2.068xl0"4 

6 

1.860xl0"4 

7 

1.672xl0-4 

8 

1.499xl0-4 

9 

1.339xl0"4 

10 

1.193xl0"4 

11 

1.059xl0-4 

12 

9.376xl0"5 

13 

8.266xl0"5 

14 

7.067xl0"5 

15 

6.034xl0"5 

16 

5.154xl0"5 

17 

4.402xl0-5 

18 

3.743xl0"5 

19 

3.183xl0“5 

20 

2.711xl0-5 

*> 

Off-Line 
.7614  iiia 

3.569xl0“4 
3 . 236xl0~4 
2 ,924xl0~4 
2.651xl0~4 
2 . 397xl0~4 
2 . 161xl0“4 
1.944xl0“4 
1.747xl0“4 
1 . 566xl0“4 
1 . 399xl0-4 
1.246xl0“4 
1.107xl0“4 
9.795xl0“5 
8 . 636xl0-5 
7 . 383x10"'’ 
6 . 305xl0-6 
5 . 385xl0"5 
4 . eooxio-5 
3.910xl0“5 
3 . 326x10"’’ 
2 . 833xl0-5 


where  Pq  =  initial  pressure 

P  =  pressure  at  another  height 
z  =  height  (km) 

H  =  scale  height  (km). 


The  scale  height  H  is  calculated  from 


ft  T 

H  =  =  6.9  km 


where  H  =  universal  gas  constant 

g  =  acceleration  of  gravity. 


Pressure  can  be  defined  through  the  ideal  gas  law  as 


P.  =  N.  k  T 
i  1  o 


where  P^^  =  pressure  at  some  altitude  i 
=  number  density  at  altitude  i 
k  =  Boltzmann  constant 


Tq  =  temperature  at  altitude  i. 


So  by  substitution  into  Equation  59 ! 


N.  k  T  =  N  k  T  exp  ( — ) 
i  o  o  o  r  H 


TT  =  exP  {hs>  •  (62) 

o 

Thus,  number  density  is  also  exponentially  distributed  with  height  in  an 
isothermal  atmosphere. 

The  lidar  equation  (see  Equation  18)  defines  atmospheric  attenuation 

due  to  gases  0  as 
g 


1  v  >  * .  ** ,  •*.  \  ■* .  *.  •*  «■,  ,■  >■  v  . ■ '  - 
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i-  /V  V  I'yy.i'. 


oig(z)  dz 


(63) 


J 


where  ZL  =  altitude  of  shuttle  lidar  (km) 

Z-r  =  height  of  thunderstorm  top  (km) 

oig  =  atmospheric  extinction  coefficient  due  to  gases  (kra-^-)  . 

It  will  be  assumed  that  c»g  follows  the  exponential  distribution  described 
by  an  isothermal  atmosphere.  Integrating  the  right-hand  side  of  Equa¬ 
tion  59,  Qg  can  be  found  by 

oZl 

Qg  W  J  ^  exp  (-  -—)  dz  (64) 


~  2T' 

.9  ’  J 


Qg  =  oig (Zt)  6.9  exp  (  - 


Consider,  ZL  =  200  km  and  Z-j  =  10  km.  Then  Zj,  -  Z-j  =  190  and  the  exponent 
of  Equation  65  approaches  one.  The  extreme  distances  separating  the 
shuttle  from  a  cloud  top  implies  that  Equation  65  can  be  rewritten  as 


Qg  ”  otg(Z<p)  6.9 


Extensive  tables  of  ag(z)  at  many  laser  wavelengths  have  been  compiled 
by  Elterman  (1968) .  In  this  thesis,  the  values  of  ag(Z)  listed  at  .8  Um 
are  used  for  both  the  on-line  and  off-line  wavelengths  (see  Table  II) . 


Attenuation  by  Aerosols 

Aerosol  particles  in  the  atmosphere  vary  in  size  from  clusters  of 
a  few  molecules  up  to  particles  of  100  ym  in  radius.  Aerosols  can  be 


be  divided  into  three  size  categories  (Pruppacher  and  Xlett,  1980) : 


H 


1.  Aitken  particles  -  particles  with  dry  radii  riO.l  um  (individual 


atoms  and  molecules) 

2.  large  particles  -  particles  with  dry  radii  0.1  <  r  <  1.0  un 
(smoke,  haze) 

3.  Giant  particles  -  particles  with  dry  radii  1.0  p  <  r  <  10  pm 
(dust) . 

This  thesis  distinguishes  between  aerosol  particles  and  cloud  parti¬ 
cles.  Aerosols  will  be  considered  as  particles  with  radii  r  <  10  pm. 

Aerosol  particles  are  injected  into  the  atmosphere  by  natural  and 
by  anthropogenic  sources.  Most  particles  originate  from  the  earth's 
surface,  interiors  of  volcanoes,  or  outer  space.  Particles  of  terrestrial 
origin  are  formed  by  gas-to-particle  conversions  or  by  mechanical  and 
chemical  disintegration  of  the  earth's  solid  and  liquid  surfaces.  Some 
typical  particles  are  made  of  soil  and  rock  debris,  sea  salt,  volcanic  debris 
sulfates,  nitrates,  and  ammonium  salts. 

Using  exponentially  derived  data,  aerosol  particle  size  distribution 
can  be  expressed  according  to  Junge's  (1972)  power  law  as 

dN(r)  =  rr(D+l) 

“dr  (67) 

where  dN(r)  =  number  density  of  particles  between  particle  radius  r  and 
dr 

C  =  normalizing  constant 

D  =  shaping  constant,  2  <_  D  <_  4,  for  particle  radii  .1  to 
10  pm. 

Elterman  (1968)  expanded  the  Junge  distribution  and  computed  extensive 

tables  of  atmospheric  extinction  due  to  aerosols  a A  at  different  wavelengths 
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from  altitudes  of  0-50  km.  His  research  examined  over  79  density 


profiles  demonstrating  how  the  aerosol  extinction  coefficient  changes 
with  height. 

Detailed  studies  by  Junge  (1972)  in  air  over  Germany  found  that  the 
total  aerosol  concentration  decreased  exponentially  up  to  about  the 
tropopause.  Above  the  tropopause,  the  aerosol  concentration  decreases 
rapidly  toward  zero  at  50  km.  The  data  supplied  by  Elterman  (1968)  at 
,8um  indicates  that  a^  is  distributed  according  to  a  power  law  from  near 
12  km  to  50  km.  A  least  squares  curve  fit  technigue  was  applied  to  this 
data  to  model  the  power  law. 

The  power  curve  can  be  expressed  as 


Z 


(68) 


where  Z  =  height  (km) 

a  =  aerosol  extinction  coefficient  (km  ^) 

A 

a,b  =  regression  coefficients. 


The  regression  coefficients  a  and  b  are  found  by  solving  the  following 


system  of  linear  equations 

( Ketter 

and 

Prawell 

n  la.  . 

Ai 

fa.,  fa2.. 

Ai  Ai 

■—  « 

a 
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•  - 
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- 

(69) 


where  n  =  number  of  data  points 


-1, 


a =  aerosol  extinction  coefficient  at  a  given  height  (km  ). 


The  coefficient  of  determination  can  also  be  expressed  in  a  numerical 


manner  by 


(70) 


2  aEZi  +  bEotAi  Zi  -  —  (XZi) 

R  =  _ n _ 

S(Zi2)  -  -  (Ezi) 2 
n 


The  model  results  are 
R2  =  .80 
a  =  2.1197 
b  =  -.25. 


Equation  68  can  now  be  rewritten  as 


or 


-.25 

Z  =  2. 1197 (aA) 


aA  =  ( 


2.1197 


) 


(71) 


(72) 


Equation  72  indicates  that  aA  (from  approximately  10  -  50  km)  decreases 
rapidly  to  almost  zero. 

Atmospheric  attenuation  due  to  aerosols  QA  is  defined  in  Equation  18 
as 

Qa  =J  aA<z)dz 

Z<p 


where  symbols  are  as  previously  defined. 

The  correct  functional  relationship  needed  to  evaluate  Qa  is  found 
by  integrating  the  right-hand  side  of  Equation  72  such  that 


“4 

( — ? - ) 

V2.1197; 


dz 


(73) 


2.1197 


Qa  =  6.73  (ZL“3  -  2t“3)  .  (74) 

Aerosol  Backscatter  Coefficient 

The  aerosol  particles  considered  in  the  above  derivation  range  from 
.1  to  10  vim  in  radius.  This  implies  that  Mie  scattering  theory  applies, 
since  the  laser  wavelength  is  on  the  same  scale  as  the  particles.  The 
aerosol  backscatter  coefficient  BA  (required  in  the  lidar  equation)  must 
be  found  using  Mie  theory  (Hinkley,  1976) . 

Investigators  typically  avoid  solving  complex  differential  equations 
in  computing  Ba.  In  general  cases,  a  simplifying  relationship  is  estab¬ 
lished  between  aA  and  BA  (Hobbs  and  Deepak,  1981) .  This  relation 
is  referred  to  as  the  backscatter-to-extinction  ratio  k  and  is  experi¬ 
mentally  measured.  The  backscatter-to-extinction  ratio  k  for  aerosols 
of  the  size  considered  in  this  thesis  was  found  by  Fernald  (1972)  to  be 

k  =  =  .03  sr_1  .  (75) 

aA 

It  is  then  a  simple  manner  to  compute  BA  from 

.03  afl  =  Ba.  (76) 


Attenuation  by  Thunderstorm  Tops 

The  top  0.5  km  of  a  thunderstorm  in  this  thesis  will  be  modeled  as 
a  cirrus  ice  cloud.  This  assumption  is  consistent  with  Byers  and 
Braham's  (1949)  model  of  a  thunderstorm  calling  for  the  upper  regions  of 


the  cumulon imbus  cloud  to  be  made  of  cirrus  type  ice  crystals.  Cirrus 


clouds  are  composed  of  several  types  of  ice  crystals  such  as  bullet-shaped 

crystals,  single  rosettes,  hexagonal  columns,  and  irregular  crystals 

(Pruppacher  and  Klett,  1980),  Ice  crystal  concentrations  range  between 
4  -3 

1  to  5x10  cm  , 

It  may  be  wondered  if  a  laser  beam  can  even  penetrate  a  cloud  without 
suffering  severe  attenuation.  Platt  (1979)  has  conducted  extensive 
research  with  lidars  probing  clouds.  He  reports  that  for  low-level  water 
clouds,  laser  penetration  is  about  100  mj  however,  for  cirrus  clouds, 
laser  penetration  is  found  to  be  >4  km.  This  is  true  since  cirrus  clouds 
have  very  low  water  liquid  water  content  this  low  internal  scattering 
from  the  water  molecule. 

Spinhirne  et  al  (1982)  used  a  high-altitude  research  aircraft  with 
a  Nd:YAG  lidar  mounted  in  the  plane  to  actually  measure  the  volume  back- 

scatter  coefficient  by  over-flying  cirrus  clouds.  The  measured  value  of 

-1  -1 

the  cirrus  cloud  top  volume  backscatter  coefficient  B  is  0.02  km  sr  . 

c 

The  extinction  coefficient  due  to  the  presence  of  cloud  particles 


ac  can  be  found  through  a  backscatter-to-extinction  ratio  k.  Platt 

(1973)  probed  cirrus  clouds  with  a  lidar  and  had  great  success  with  a 

value  k  =  .05  sr  Accordingly,  this  thesis  will  use  the  following 

values  for  B  and  a  : 

c  c 

B  =  .02  km  1  sr  1 
c 

Bo  .02  km"1  sr-1  .  .  -1 


The  only  other  assumption  that  investigators  make  is  that  Bc  and 
ac  will  be  constant  throughout  the  depth  of  cloud  being  probed  (1/2  km 


in  this  thesis) . 

Atmospheric  attenuation  due  to  cloud  particles  Qc  is  defined  in 
Equation  18  as 

p  ZT 

Qc  =  \  Oj,  dz  =  .5  ac.  (77) 

ZT  -1/2  km 

Miscellaneous  Cloud  Top  Parameters 

With  knowledge  of  ac,  it  is  possible  to  compute  other  cloud  top 
parameters  such  as  cloud  particle  number  density  Nc  and  the  liquid  water 
content  Wc  existing  on  ice  crystals.  The  number  density  is 

Nc  =  Bc  (78) 

27Ta^ 

where  a  =  cloud  particle  radius  (a  typical  value  is  50  pm) . 

The  cloud  water  content  can  be  found  by 

Nc  =  -|  it  a2  p  Nc  (79) 

where  p  =  cloud  particle  mass  density  (.92  gm  cm-2  for  ice  clouds, 
Russell,  1982) . 

Attenuation  by  Ozone 

The  extinction  due  to  ozone  a 3  is  a  strong  function  of  altitude 
z  and  wavelength  X.  Elterman  (1968)  defines  ou  as 


a3  =  Ay  (X)  d3  (2) 


(80) 


where  Av  =  pure  ozone  absorption  coefficient  (cm-*) 

D3  =  ozone  equivalent  thickness  (cm  km"1) . 

Av  for  the  wavelength  used  in  this  thesis  is  l.OOxlO-^  cm-*. 

Values  for  D3  have  been  tabulated  by  Elterman  (1968)  for  heights  of  0-50 
km. 

This  thesis  is  concerned  with  ozone  03  at  heights  from  20-50  km. 
Elterman  (1968)  shows  that  a3  (20  km,  .8  pm)  =  1.64xl0-4  km-1  while 
a3  (50  km,  .8  ym)  =  1.86x10“^  km-*,  since  these  values  are  extremely 
small  compared  to  oig,  aA,  and  ac  the  effects  of  ozone  are  not  considered 
in  this  thesis.  Figure  8  is  a  profile  of  a3  at  .8  in  as  a  function  of 
height  (Elterman,  1968). 


-»■ .  rw'  tr  <  <r !  in  v  v  rs-vw  ■vjirjT'JF'rw'; 


VI.  Laser  Considerations 

Laser  Beam  Divergence 

As  a  laser  beam  leaves  the  optical  cavity  and  propagates  it  will 
slowly  diverge  because  of  diffraction.  Assuming  the  laser  beam  is  a 
Hermit-Gaussian  type,  beam  divergence  0  can  be  expressed  as 

0  =  .J-i— 

IT  W0 

_  4  X 
IT  2W0 

1.27  (81) 

2W0 

where  0  =  full  angle  beam  divergence 

X  =  laser  operating  wavelength 

WQ  =  beam  spot  size  (actually  this  is  the  beam  radius  value) 
at  beam  waist  (O'Shea  ot  ad,  1978) . 

Another  factor  of  interest  to  the  laser  user  is  spot  size  on  the 
target  W.  This  thesis  is  interested  in  aiming  the  laser  beam  at 
the  rising  towers  protruding  above  the  anvil  floor  so  it  would  be 
desirable  for  W  to  be  less  than  the  diameter  of  an  updraft  tower. 
Typical  thunderstorm  updraft  radius  range  from  1-3  km  (Kropfli  and 
Miller,  1976)  . 

A  Hermit-Gaussian  beam  with  a  spot  size  W  on  a  target  at  a 
distance  z  can  be  expressed  according  to  Nielsen  (1980)  as 


where  z  is  the  distance  from  the  minimum  spot  size  (essentially  the 


laser's  front  end)  to  the  target.  It  must  be  noted  that  Equation  82 
is  valid  only  when  the  target  which  in  this  case  is  a  storm  top  lies 
in  the  laser's  far  field.  A  laser's  near  field  is  that  distance  z 
where  the  beam  remains  well  collimated.  The  Rayleigh  range  can  now 
be  defined  as  that  distance  z  where  the  laser  beam  begins  to  spread. 
The  Rayleigh  range  Zr  in  equation  form  is 


(83) 


Once  the  target  distance  Zi>  exceeds  20-100  times  Zr,  the  laser  beam 
enters  the  far  field  region  and  Equation  82  becomes  effective. 

Figure  9  illustrates  the  situation  considered  in  this  thesis. 

It  will  now  be  demonstrated  that  a  storm  top  lies  in  the  laser's 
far  field  and  W  at  the  storm  top  will  be  found.  This  example  will 
assume  that  the  lidar  is  in  a  200  km  orbit  with  a  thunderstorm  top 
at  10  km.  Spot  size  on  target  can  be  found  using  Equation  82  re¬ 
written  here  for  convenience 


W(Z)  =  W0 


1/2 


W0  =  laser  beam  radius  at  instrument  aperture  (assume  a  typical  value 
of  1  cm  (Nielsen,  1980) 

W(Z)  3  beam  radius  at  target 

A  =  laser  operating  wavelength  ( .7696x10"^  m) 

z  *  distance  from  lidar  to  storm  top  (190  km  here) . 


To  use  this  formula,  the  target  must  lie  in  the  far  field.  The  far 
field  check  is  performed  by  testing  if  z  is  20  to  100  times  greater 


than  Equation  83: 


60 


far  field  check 


z  >  100  times 


X 


_ _  ^  tt(.01  m)2 

190,000  m  .7696xl0“6m 

190,000  m  >  (100)  (408  m) 
190,000  m  >  40,800  m. 


So  target  (storm  top)  is  in  the  far  field. 

The  beam  radius  on  target  W(190  km)  according  to  Equation  82 
is  4.655  m.  The  circular  area  illuminated  at  the  cloud  top  is 
68.059  m2. 

These  values  are  much  smaller  than  the  average  updraft.  Such 
values  indicate,  for  the  first  time,  an  orbiting  sensor  could  observe 
cloud  top  structure  on  a  scale  previously  obtainable  only  by  high 
altitude  aircraft. 


Thermal  Blooming 

The  atmosphere  absorbs  photons  from  the  beam  as  it  propagates. 
This  absorbed  energy  is  not  lost  but  instead  heats  and  expands  the 
channel  of  air  surrounding  the  beam.  This  rise  in  temperature  of  the 
air  channel  is  thermal  blooming,  and  the  effect  is  to  cause  the  beam 
to  spread  farther  than  Equation  82  predicts. 

Thermal  Blooming  is  greatly  reduced  by  pulsing  the  beam,  since 
a  short  pulse  does  not  stay  in  one  location  long  enough  for  the  air 
channel  to  expand.  The  time  scale  of  thermal  blooming  can  be  found 


according  to  Nielsen  (1980)  as 


where  t_.  =  time  required  for  thermal  blooming  to  occur 
W  =  beam  radius  at  specified  distance 
j  =  speed  of  sound  (3x10^  cm/sec) . 

The  previous  example  found  W(190)  =  4.66  m.  This  implies  tj  =  1.55x 
10-2  sec.  The  pulse  width  used  in  this  thesis  is  3.33xlO~®  sec.  So 
thermal  blooming  does  not  have  time  to  occur  and  will  therefore  be 
neglected. 

Dye  Lasers 

Dye  lasers  must  be  used  when  DIAL  experiments  are  being  per¬ 
formed  such  as  the  DIAL  procedure  employed  by  this  thesis.  DIAL 
experiments  require  two  wavelengths:  one  at  the  absorption  line 
of  the  molecule  studied  and  a  second  nearby  wavelength  where  ab¬ 
sorption  by  the  molecule  is  negligible.  Furthermore,  the  absorption 
lines  of  most  molecules  are  not  at  standard  laser  wavelengths. 

Dye  lasers  can  meet  all  DIAL  experimental  requirements  since  they 
can  be  tuned  to  non-standard  wavelengths.  Dye  lasers  can  be  frequency 
tuned  to  any  frequency  in  the  visible  spectrum.  Such  lasers  work 
well  in  the  pulse  mode  and  are  characterized  by  high  output  power, 
high  gain,  and  low  cost. 

The  active  medium  of  dye  lasers  consist  of  an  organic  dye 
dissolved  in  a  solvent.  Industry  is  developing  a  dye  laser  that 
does  not  require  a  dye  cell.  Instead  a  rotating  solid  plastic  disc 
is  impregnated  with  dye  (Greco,  1980) .  A  suitable  dye  that  can  be 
used  to  reach  the  wavelengths  chosen  in  this  thesis  (both  wavelengths 
centered  near  .76  yim)  is  nilebluo--A  perchlorate  (O’Shea  et  al. 


1978)  or  oxazine  725  dye  in  ethanol  (Greco,  1980) . 

Dye  lasers  are  optically  pumped  which  implies  that  another  laser 
can  be  used  to  pump  the  dye  laser.  The  only  requirement  is  for  the 
pumping  laser  to  have  an  output  near  the  peak  dye  absorption  band. 

A  dye  laser  operating  in  the  oxygen  "A"  band  could  be  pumped  by 
nitrogen  or  krypton  lasers.  The  krypton  laser  is  chosen  as  the 
pumping  laser  in  this  thesis  since  Kalshoven  et  jLL  (1980)  had  great 
success  with  such  a  system  while  operating  in  the  oxygen  "A”  band. 


VII.  Instrumental  Considerations 

Integration  Time 

The  integration  time  tj  will  be  defined  as  that  time  interval 
for  which  the  lidar  receiving  telescope  is  turned  on  or  receiving 
signal  and  noise  photons.  This  integration  time  for  Kts  pulses 
can  be  expressed  as 


tl  = 


2dKts 

c 


(85) 


where  d  =  vertical  resolution  (km) 

K  =  number  of  pulses  per  second  (PPS) 
ts  =  sounding  time  (see  below) 
c  =  speed  of  light  (km  sec-^-)  . 


The  integration  time  can  be  calculated  using  the  correct  values 
from  the  list  of  laser  characteristics  found  in  Table  V. 


Sounding  Time 

The  sounding  time  ts  is  defined  as  that  time  interval  required 
to  place  1000  shots  (Kt^  pulses)  into  a  cloud  top.  Sounding  time 
will  be  taken  as  20  seconds  (20  sec  x  50  shots  sec”^-  =  1000  shots)  . 


Laser  Power 

Huf faker  (1978)  indicates  that  the  space  shuttle  can  allocate 
a  maximum  of  3000  watts  to  a  shuttle-borne  lidar.  The  lidar  system 
developed  for  this  thesis  is  a  three  laser  system  where  two  of 


these  are  dye  lasers.  The  first  dye  laser  is  tuned  to  an  02  ab- 
sorption  line,  the  second  dye  laser  is  tuned  to  a  nearby  spectral 


Table  V 


Laser  Characteristics 


Energy  per  Pulse  (Ep) 

6J 

Laser  Output  power 

300  watt 

Peak  Power  (Pp) 

1.8  MWatt 

‘Shuttle  Power  available  to  Lidar 

3000  Watt 

Lidar  System  Efficiency 

10% 

Pulse  Width  (tp) 

3 .33x10 “6 

Vertical  Resolution  (d) 

0.5  km 

‘Receiver  Diameter 

1.25  m 

On-line  Wavelength 

.7696  ym 

Off-Line  Wavelength 

.7614  ym 

‘Optical  Filter  Bandpass  (AX) 

1  nm 

‘Receiver  FOV 

1  mrad 

Pulse  Repetition  Rate  (k) 

50  pps 

Integration  Time  tj 

3.33xl0-3 

Sounding  Time  ts 

20  sec 

indicates  values  recommended  by  Greco  (1980) 
indicates  value  recommended  by  Huffaker  (1078) 


position  where  1>2  absorption  is  negligible.  The  third  laser  (a 
krypton  gas  laser)  will  optically  pump  the  first  two.  Such  a  com¬ 
plex  system  of  three  lasers  along  with  a  well  designed  optical  sys¬ 
tem  is  generally  considered  to  be  about  10%  efficient  (Huf faker, 
1978)  .  Therefore,  laser  output  power  could  be  as  large  as  300 
watts . 

The  energy  per  pulse  using  this  power  level  is  calculated  as 


EP  "  K 


(86) 


where  Ep  =  energy  per  pulse  (J) 

PA  =  power  available  to  laser  (watts) 
e  =  laser  system  efficiency (see  Eg.  14) 


Using  Equation  86,  Ep  will  be  taken  as  6j  per  pulse.  Commercial 
lasers  developing  6j  per  pulse  will  be  available  by  the  mid  eighties 
according  to  Huffaker  (1978)  and  Greco  (1980) . 

Peak  Power 

High  powered  pulsed  lidars  operated  from  the  shuttle  will  de¬ 
velop  very  high  peak  power.  High  peak  powers  could  represent  a 
danger  to  the  human  eye  for  an  observer  viewing  the  laser  directly 
from  the  earth's  surface.  Eye  safety  is  discussed  in  Chapter  X  of 
this  thesis .  Peak  power  Pp  can  be  found  by 


where  tp  =  the  pulse  width  (see  Table  V) . 
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Vertical  Resolution 


The  vertical  resolution,  or  depth  of  the  volume  sensed,  will 
be  taken  as  0.5  km.  Penetrating  a  thunderstorm  top  beyond  0.5  km 
is  thought  to  be  impossible  since  the  laser  beam  would  be  greatly 
attenuated  by  the  presence  of  super-cooled  water  droplets  suspended 
in  the  updraft. 

Another  problem  can  occur  if  the  beam  could  penetrate  the  cloud 
top  much  beyond  0.5  km.  The  laser  beam  will  encounter  warmer 
(243°K)  temperatures  as  penetration  increases.  Cloud  microphysics 
predicts  hexagonal  ice  plates  in  clouds  at  a  temperature  range  of 
243-261°K  (Pruppacher  and  Klett,  1980) . 

Platt  (1978  a,b)  showed  that  super,  large  backscatter  returns 
are  possible  if  lidar  signals  are  reflected  off  horizontally 
oriented  hexagonal  ice  plates.  This  phenomenon  occurs  only  if  the 
zenith  angle  ip  is  zero.  Figure  10  illustrates  this  situation  with 
a  crystal  of  thickness  t.  To  make  manners  worse,  these  plate 
shaped  crystals  fall  naturally  with  their  long  axis  horizontal. 

Platt  (1978  a)  further  explains  that  sudden  high  increases  in 
the  signal-to-noise  ratio  (see  this  chapter)  are  possible  when 
laser  light  is  reflected  off  even  a  single  horizontally  oriented 
hexagonal  ice  plate  as  the  lidar  is  scanned  across  the  vertical 
(zero  zenith  angle) .  Therefore,  it  will  be  assumed  that  no  such 
ice  crystals  exist  in  the  colder  top  0.5  km  of  a  thunderstorm. 

A  final  barrier  that  could  stop  a  laser  beam  from  ever  reach¬ 
ing  the  storm  top  is  for  dense  cirrus  clouds  to  hide  the  thunder¬ 
storm  altogether.  Such  cirrus  are  produced  by  previous  convection. 


This  canopy  of  cirrus  over  thunderstorms  does  occasionally  form 


(Adler  and  Fenn,  1979c) .  This  problem  is  disregarded  in  this  thesis 
The  vertical  resolution  d  can  be  computed  by 

d  =  2_|E  (88) 

where  symbols  are  as  previously  defined. 

Background 

The  onboard  receiving  telescope  will  receive  return  rad  Lation 
from  thunderstorm  ops  and  from  the  background.  Three  observations 
will  be  considered: 


a.  sunlit  cloud  with  albedo  A  =  .8 

b.  sunlit  ocean  with  albedo  A  =  .1 

c.  moonlit  cloud  with  albedo  A  =  .8. 


The  albedo  A  is  the  percentage  of  radiation  reflected  by  a  body. 

Background  radiance  (watt  m“2  st-^  ym“i)  can  be  found  ac¬ 
cording  to  Russell  (1982)  as 


LX  = 


AE\  (X)  cos 


where  E\  (X)  =  solar  spectral  irradiance  (watt  m~2  ym-^) 
\p  =  zenith  angle. 


The  zenith  angle  ip  is  chosen  to  be  zero  so  that  only  worse  case 
situations  are  considered  such  as  local  noon  with  the  sun  directly 
overhead  or  the  full  moon  directly  overhead  at  night.  The  value  of 
E\  (X)  at  .76  ym  will  be  1211  watts  m-2  ym“l  (Thekaekara,  1974) . 


Table  VI  displays  upward  spectral  radiance  for  a  downward  pointing 
lidar  operating  near  .76  um.  Values  in  Table  VI  reflect  worse  case 
background  situations.  Moonlit  values  correspond  to  the  ratio: 


E^.  (sun>  =  10-6 
Ex  (moon) 


according  to  Russell  (1982) . 


(90) 


Table  VI 

Upward  Spectral  Radiance  (watts  m~2  nm”1  sr“l)  at  .76  pm 


Sunlit 

Moonlit 

Cloud  top 

3.084xl0-1 

3.084x10 

Ocean 

3.855xlO~2 

3.855x10 

Signal-to-Noise  Ratio 

The  power  received  by  the  lidar' s  receiver  PL  given  by  Equa¬ 
tion  18  contains  three  sources  of  noise  which  collectively  limit 
the  sensitivity  of  the  receiver.  Noise  results  from  the  random 
arrival  of  photons  from  cloud  induced  laser  returns  and  the  back¬ 
ground.  This  random  arrival  is  described  by  a  Poisson  distribu¬ 
tion.  Another  source  of  noise  is  the  random  generation  of  dark 
current.  This  dark  current  occurs  within  the  photomultiplier  even 
if  no  photons  are  incident.  It  can  be  neglected  in  photomultipliers 
if  they  are  cooled  to  liquid  nitrogen  temperatures  (77°k) . 

Investigators  typically  discuss  noise  in  terms  of  photo¬ 
electrons  produced  by  the  detector.  Accordingly,  the  signal-to- 
noise  ratio  S/N  can  be  written  as  (Smith  and  Platt,  1978)  : 


(n  +  n  ) 
s  b 


where  =  signal  count  from  cloud  top 


n,  =  background  count, 
b 


Signal  counts  are  calculated  by 


n  =  7“  e  P  t 
s  he  LI 


where  X  =  laser  wavelength  (m) 

-34 

h  =  Planck's  constant  (6,63x10  J-sec) 

8  -1 

c  =  speed  of  light  (3.0x10  m  sec  ) 
e  =  system  efficiency  (.10) 

P  =  received  power  collected  by  receiver  optics  (watts) 
Li 

tj  =  integration  time  (3.33x10  3  sec). 


Included  in  e  is  the  optical  system  efficiency  and  the  receiving 

he 

system  (photomultiplier)  efficiency.  The  quantity  —  is  the  energy 
per  photon . 

Background  counts  are  calculated  by 


nb  = 


— —  e  P,  t 
he  b  I 


where  Pfa  (watts)  is  the  power  from  the  cloud  background  and  is 
calculated  by 


P,  =  fiAX  A  L 
b  r 


where  receiver's  field  of  view  (sr) 


AX  =  receiver's  optica)  bandpas 


f i Iter  (nm) 


VIII.  System  Results 


Scenario  I 

The  lidar  equation  (see  Equation  18)  can  now  be  solved  using  the 
quantities  defined  in  the  atmospheric  considerations  chapter  and  the 

I 

laser  system  characteristics  specified  in  Table  V. 

Tables  VII  -  XIII  display  the  numerical  results  of  solving  the 
lidar  equation  for  varying  shuttle  altitudes  and  thunderstorm  heights 
Signal-to-noise  calculations  are  performed  for  situations  of  a  worse 
case  daylight  (local  noon)  and  a  worse  case  at  night  with  a  full 
moon  overhead. 

The  results  shown  in  Tables  VII  -  XIII  indicate  that  the  shuttle- 
borne  lidar  proposed  in  this  thesis  can  successfully  probe  the  top 
0.5  km  of  a  thunderstorm  and  return  a  measurable  signal.  The  night¬ 
time  values  of  S/N  are  very  good  as  would  be  expected  with  reduced 
background  radiance  from  the  storm  tops. 

Two  scenarios  will  be  considered  in  this  thesis.  This  first 
scenario  will  demonstrate  the  shuttle  lidar' s  ability  to  determine  a 
single  temperature.  At  the  end  of  the  following  chapter  on  thunder¬ 
storm  vertical  velocities,  a  second  scenario  will  consider  a  second 
temperature  measurement  of  the  same  area  at  a  later  time. 

Consider  a  scenario  where  the  shuttle  is  orbiting  at  200  km. 
Suppose  that  through  coordination  with  meteorologists  on  the  ground 
using  weather  radar,  the  shuttle  lidar  is  aimed  at  an  overshooting 
top  at  10  km  above  the  surface.  Using  the  results  shown  in  Table  XI, 
Equation  19  can  be  solved  for  the  resonant  absorption  coefficient 
dg0  such  that 
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_ (1.439)  (1421.436  cm  ~) _ 

Ti+1  In  £(  1742.73°KJ/'i  itfl)  (.9999)  (235°K-3/2) 

1  (10)  { .7486)'  .0512  km  x 

=  235.09°K. 

If  needed,  further  iterations  could  be  performed  using  the  previous 
result  for  as  the  seed  temperature  until  successive  results  give 
the  same  results  to  within  0.1°K, 

Attitude  Control 

The  space  shuttle's  ability  to  hold  a  navigation  axes  oriented 
to  nadir  is  0.5  degrees  (3  sigma)  up  to  one  hour  (Huffaker,  1978). 

In  addition,  there  is  another  two  degrees  of  error  possible  due  to 
misalignment  of  the  navigation  subsystem  to  the  cargo  bay.  Therefore, 
shuttle  aiming  accuracy  for  a  lidar  mounted  in  the  cargo  bay  is  -  2.5 
degrees. 


Data  Collection  Intervals 

The  space  shuttle  operates  in  low  earth  orbit  typically  at  200  km 
with  an  orbital  velocity  of  5  km  sec  This  implies  that  the  time 
interval  during  which  a  storm  top  remains  within  range  of  the  lidar 
is  limited.  Indeed,  this  short  loitering  time  is  a  major  factor  for 
anv  orbit  except  a  geosynchronous  one. 

It  must  be  determined  how  many  data  collection  intervals  or  how 
many  20  second  sounding  periods  are  available  to  the  instrument  during 
a  typical  flyover.  Figure  11  illustrates  a  flyover  with  a  thunderstorm 
top  at  12  km.  It  can  be  seen  that  about  four  20  second  data  collection 
intervals  are  possible  before  the  instrument  becomes  sicxnal-to-noise 


Table  VII 


Tabulated  Results  from  Solving  Lidar  Equation 

Shuttle  Altitude  =  100  km 


Storm  Top 
(ft) 

Return 

Power 

On-Line 

(nw) 

Return 

Power 

Off-Line 

(mw) 

S/N 

Daylight 

S/N 

Moonlight 

32810 

1.497xl0-4 

9.86 

438.78 

36091 

1. 808xl0-4 

3 .  HOxlO-4 

11.90 

482.23 

1^1 

2 . 171xl0-4 

3 . 213xl0-4 

14.29 

528.51 

42653 

2 . 466xl0-4 

3 . 260xl0-4 

16.24 

563.39 

45934 

2 . 661xl0-4 

3 . 338xl0-4 

17.52 

585.25 

49215 

2.904xl0“4 

3.418xl0~4 

19.11 

611.34 

52496 

3 . 195xl0-4 

3 . 500xl0-4 

21.03 

641.32 

55777 

3 . 561xl0-4 

3 . 586xl0-4 

23.44 

677.08 

Stored  Top 
(km) 

Storm  Top 
(ft) 

Return 

Power 

On-Line 

(yw) 

Return 

Power 

Off-Line 

(MW) 

S/N 

Daylight 

S/N 

Moonlight 

10 

32810 

9.170xl0-5 

1 . 862xl0-4 

6.04 

343.18 

11 

36091 

1 . 102xl0-4 

1.896xl0-4 

7.25 

376.28 

12 

39372 

1 . 317xl0-4 

1 . 949xl0-4 

8.67 

411.40 

13 

42653 

1.488xl0~4 

1. 967xl0-4 

9.80 

437.46 

14 

45934 

1 . 598xl0~4 

2 . 004xl0-4 

10.52 

453.27 

15 

49215 

1 . 734xl0~4 

2.041xl0"4 

11.41 

472.24 

16 

52495 

1 . 898xl0~4 

2 . 079xl0-4 

12.49 

494.07 

17 

55777 

2 . 118xl0-4 

13.85 

520.20 

— 

Return  Return 

Power  Power 

Storm  Top  Storm  Top  On-Line  Off-Line  S/N 

(km)  (ft)  (yw)  ( yw)  Daylight 

10  32810  6.187xl0~5  1.256xl0-4  4.07 

11  36091  7.411xl0-5  1 . 275xl0-4  4.88 

12  39372  8.828xl0-5  1.307xl0~4  5.81 

13  42653  9 . 947xl0~5  1.315xl0-4  6.55 

14  45934  1.064xl0“4  1.335xl0~4  7.01 

15  49215  1 . 151xl0~4  1.355xl0-1  7.58 

16  52496  1. 256xl0-4  1.375xl0"4  8.27 

17  55777  1 . 387xl0-4  1.397xl0-4  9.13 


S/N 

Moonlight 

281.68 

308.40 

336.69 

357.46 

369.78 

384.62 

401.74 

422.26 


Return 

Power 

Storm  Top  Storm  Top  On-Line 
(km)  (ft)  (yw) 

10  32810  3.359xlO-5 

11  36091  4 . 009xl0~5 

12  39372  4 . 756xl0-5 

13  42653  5 . 339xl0-5 

14  45934  5.689xl0-5 

15  49215  6.130xl0-5 

16  52496  6.659xl0~5 

17  55777  7 . 325xl0-5 


Return 

Power 

Off-Line  S/N  S/N 

(yw)  Daylight  Moonlight 


6 . 821xl0-5 

2.21 

207.13 

6.897xl0"5 

2.64 

226.43 

7.040xl0-5 

3.13 

246.79 

7 . 056xl0-5 

3.52 

261.54 

7 . 136xl0-5 

3.75 

270.05 

7 . 215xl0-5 

4.04 

280.35 

7.295xl0-5 

4.38 

292.26 

7 . 377xlO~5 

4.82 

306.60 

Stora  Top 
(ft) 

32810 

36091 

39372 

42653 

45434 

49215 

52496 

55777 


Return 

Power 

On-Line 

(uw) 

Return 

Power 

Off-Line 

(uw) 

S/N 

Daylight 

S/N 

Moonlight 

2.624xl0“5 

5 . 327xl0-5 

1.73 

182.82 

2 . 127xl0-5 

5 . 379xl0-5 

2.06 

199.77 

3.705xl0~5 

5.485xl0“5 

2.44 

217.63 

4.154xl0-5 

5.490xl0-5 

2.74 

230.52 

4 . 421xl0”5 

5 . 545xlO~5 

2.91 

237.87 

4. 757xl0-5 

5. 59Sxl0~5 

3.13 

246.80 

5.161xl0-5 

5.654xl0~5 

3.40 

257.14 

5.670xl0-5 

5 . 710x10”'’ 

3.73 

269.59 

Return  Return 

Power  Power 

Storm  Top  Storm  Top  On-Line  Off-Line  S/N  S/N 

(km)  (ft)  (uw)  (uw)  Daylight  Moonlight 

10  32810  2.105xl0~5  4.275xl0'5  1.34 _ 163.56 

11  36091  2 . 507xl0~5  4.313xl0~5  1.65  178.67 

12  39372  2 . 968x10“ 5  4.393xl0~5  1.95  194.58 

13  42653  3 . 324xl0~5  4.393xl0-5  2.19  206.02 

14  45934  3.534xl0*5  4.432xl0-5  2.33  212.50 

15  49215  3 . 799xl0~“  4.471xl0-5  2.50  220.38 

16  52496  4.117xl0"5  4.510xl0-5  2.71  229.50 

17  55777  4.519xl0~5  4.550xl0~5  2.98  240.50 


shuttle 

direction 


100  km  100  km  100  km 


shuttle  altitude 

200 

km 

shuttle  velocity 

5 

km  sec' 

thunderstorm  top 

12 

km 

sounding  time 

20 

sec 

limited  or  the  aiming  angle  Y2  becomes  too  small  for  shuttle  guidance 
systems  to  properly  aim  the  instrument.  At  night,  S/N  is  very  good 
and  more  than  4  data  collection  intervals  are  possible,  but  the  pointing 
problem  could  be  a  factor.  Figure  12  shows  the  same  viewing  scenario 
except  now  the  sounding  time  t  is  10  sec  instead  of  20  sec  as  in  Figure 
11.  There  are  more  intervals,  but  this  does  not  imply  that  several 
thunderstorm  tops  can  be  surveyed. 

The  distance  between  the  shuttle  and  a  storm  top  Z  is  found 
from  the  expression 

ZU=1ZLl+d-2ZLia“S  VV2  (95) 

where 

ZT ,  =  shuttle  altitude  above  storm  top  at  start  time  t 
Ll  o 

Z„ .  =  shuttle  altitude  above  storm  top  at  t  +  x  sec 
L3  ^  o 

d  =  distance  shuttle  traveled  at  5  km  sec*"^ 

Y  =  inclusive  angle  between  Z  .  and  d  in  degrees . 

The  time  scale  for  temperature  changes  in  storm  tops  reflecting 
changes  in  intensity  is  in  the  range  of  several  degrees  per  minute 
(Kack  et  al,  1982).  The  shuttle-borne  lidar  would  be  within  range  of 
one  thunderstorm  top  for  80  to  100  seconds  (see  Figure  11)  before  range 
reduces  S/N  or  the  aiming  angle  becomes  too  critical  to  ensure  the  same 
spot  is  being  measured.  Table  XIV  shows  values  of  S/N  for  different 
values  of  sounding  time  t  with  the  shuttle  orbiting  at  200  km.  Storm 
tops  from  10-  17  km  were  probed  with  t  =  10,  20,  and  30  seconds  for 
a  worse  case  daylight  situation. 


shuttle 

direction 


8' 


It  is  obvious  that  S/N  increases  as  t  increases  as  noted  by 

s 

Equation  85.  Since  a  minute  is  required  between  temperature  measure¬ 
ments,  t  =20  sec  will  provide  the  largest  S/N  and  ensure  the  cloud 
s 

top  stays  within  range.  Thus,  one  storm  top  can  be  measured  within  a 
loitering  time  of  100  seconds  if  a  rate  of  change  in  temperature  is 
required.  A  maximum  of  4  to  5  temperature  measurements  axe  possible 
if  only  single  measurements  are  required  within  80  to  100  seconds. 

A  possible  operating  sequence  where  t  =  20  sec  is  for  the  lidar 
to  lock  onto  a  desired  storm  top  and  determine  a  temperature  which 
would  require  20  seconds.  A  minute  later,  a  second  measurement  would 
be  taken  and  the  rate  of  change  in  temperature  computed  (Scenario  II 
in  Chapter  IX  will  do  this).  The  lidar  could  then  be  aimed  at  the  next 
thunderstorm  of  interest. 
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IX.  Thunderstorm  Vertical  Velocity 


Vertical  velocity  is  a  fundamentally  important  variable  which 
is  required  by  all  numerical  weather  prediction  models;  however, 
vertical  velocities  are  not  usually  measured  directly  and  are  there¬ 
fore  obtained  indirectly  with  empirical  relations. 

Adler  and  Fenn  (1979  c)  studied  thunderstorm  vertical  velocities 
with  5  minute  interval  infrared  geosynchronous  satellite  data.  They 
computed  vertical  velocities  by  an  adiabatic  method  expressed  as 


1  d  Tbb 
st  at 

3Z 


(96) 


where  w  =  vertical  velocity  of  cloud  top 

3T/3Z  =  a  reasonable  vertical  lapse  rate  (8°k/km) 
Tbb  =  minimum  black-body  cloud-top  temperature. 


Adler  and  Fenn  (1979  c)  demonstrated  that  the  intensity  of  a 
thunderstorm  is  correlated  with  the  occurrence  of  severe  weather  on 
the  ground.  They  discovered  that  the  first  report  of  a  tornado  took 
place  during  or  just  after  a  period  of  cloud  top  ascent.  This  was 
apparent  in  seven  out  of  eight  cases  studied  with  an  infrared  geo¬ 
synchronous  satellite  field  of  view  (FOV)  of  10  km  on  a  side.  Even 
with  this  poor  spatial  resolution,  it  was  determined  that  a  30  minute 
lead  time  could  be  provided  to  the  public  of  a  potential  tornado. 

Figure  13  illustrates  a  time  vs  temperature  trace  for  two  tornado 
production  thunderstorms.  These  storms  are  part  of  a  large  complex 
which  existed  on  May  6,  1975  and  covered  an  area  stretching  from 


South  Dakota  to  Texas. 


Storm  top  temperatures  wore  measured  by  Adler 
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and  Fenn  (1979  a)  using  an  infrared  geosynchronous  satellite.  The 


smoothed  traces  in  Figure  13  illustrate  the  theory  that  severe  weather 
occurs  during,  or  soon  after,  a  period  of  cloud  top  ascent  indicating 
the  storm  is  growing  and  cooling. 

When  studying  thunderstorm  vertical  velocities,  a  distinction 
must  be  made  between  the  vertical  velocity  of  the  cloud  top  and  the 
vertical  velocity  of  the  updraft.  Due  to  mixing  of  air  at  the  cloud 
top,  there  is  a  difference  between  the  two  velocities.  Kyle  et  ad 
(1976)  flew  aircraft  into  updrafts  and  was  able  to  fit  a  relationship 
to  the  results  which  can  be  written  as 


w  =  wQ  exp 


-a<f> 


(97) 


where  w  =  updraft  velocity  at  distance  r 
wQ  =  peak  updraft  velocity 
a  =  regression  constant  2.3 
r  =  radial  distance  from  updraft  core 

R  =  radius  cf  the  updraft  1-3  kra  (Kropfli  and  Miller,  1976) . 


Adler  and  Fenn  (1979  c)  calculated  an  average  vertical  velocity 


w  over  a  circular  area  of  radius  r  .  The  result  is 
2 


—  w0  ,R  . 
w  =  — -  (— T-) 
a  r1 


1  -  exp 


rl  2 

-  a 


(98) 


This  equation  links  the  vertical  growth  of  a  thunderstorm  to  its  up¬ 
draft  velocity. 


Adler  and  Fenn  (1979  a,b,c)  used  the  SMS-2  and  GOES-1  infrared 
satellites  which  have  FOVs  of  100  km2  (10  km  on  a  side) .  This  implies 
that  the  -..rcular  radius  of  integration  r^  is  5.6  km.  A  typical 


updraft  radius  R  is  3  km  according  to  Kropfli  and  Miller  (1976) . 

Through  Equation  98,  it  is  seen  that 

w  =  .12  wQ. 

This  implies  that  the  updraft  velocity  of  a  thunderstorm  is  8.3  times 
larger  than  the  storm's  mean  vertical  velocity  as  seen  by  a  satellite. 
It  is  now  quite  evident  that  vertical  velocity  measurements  are  scale 
dependent. 

Earlier,  it  was  determined  that  a  shuttle-borne  lidar  could  il¬ 
luminate  68.059  m^  of  cloud  top.  Using  R  =  3  km  and  r^  =  4.655  m  in 
Equation  98,  it  is  seen  that 

w  &  1  wQ. 

Vertical  velocities  can  now  be  measured  on  the  same  scale  as  the 
actual  updraft.  For  the  first  time,  individual  thunderstorm  updrafts 
can  be  remotely  sensed  from  orbit. 

It  is  important  to  note  that  the  lidar  is  actually  measuring  cloud 
top  temperature  and  only  through  Equation  96  can  the  vertical  velocity 
be  determined. 

Equation  96  can  be  rewritten  for  updraft  velocities  measured  by 
the  lidar  as 

1  dT 

W  =  8  ‘ 3  Jrr  dF  <"> 

9Z 

Thunderstorm  top  vertical  velocities  reported  by  Adler  and  Fenn 
(1979  c) ,  for  example,  can  now  be  altered  to  appear  as  if  they  were 
lidar  derived  (see  Figures  14  and  15). 


vertical  velocity  m  sec 


Figure  14  is  a  composite  profile  of  severe  and  non-severe  thunder¬ 


storms  observed  in  a  storm  complex  stretching  from  Oklahoma  to  Missouri 
on  24  April  1975.  The  numbers  in  parentheses  are  the  number  of  cases. 
Both  severe  and  non-sevcre  storms  slow  in  their  ascent  as  they  approach 
the  tropopause. 

Figure  15  is  similar  except  the  storm  complex  runs  from  Texas  to 
South  Dakota  and  occurred  on  6  May  1975.  Table  XV  illustrates  some 
statistical  differences  between  severe  and  non-severe  thunderstorms 
for  the  two  storm  complexes. 


Table  XV 

Statistics  on  24  April  and  6  May  1975  Storm  Complexes 


24  April  1975  Storm  Complex 


Severe  Storms : 

Average  Vertical  Velocity 

28.86  msec"1 

@  229°k 

Standard  Deviation 

12.06  msec"1 

Non-severe  Storms: 

Average  Vertical  Velocity 
Standard  Deviation 


15.54  msec"1  @  231°k 
6.81  msec"1 


6  .May  1975  Storm  Complex 

Severe  Storms: 

Average  Vertical  Velocity  24.59  msec"1  @  223°k 

Standard  Deviation  9.80  msec"1 


Non-severe  Storms: 

Average  Vertical  Velocity 
Standard  Deviation 


15.62  msec"1  0  227°k 
7.53  msec  1 
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The  t-distribution  test  can  be  performed  to  see  if  there  is  a 


difference  between  vertical  velocity  of  severe  and  non-severe  storms 
for  each  storm  complex. 

The  independent-t-test  assumes  the  mean  vertical  velocity  of 
severe  and  non-severe  storms  are  normally  distributed  with  variances 
unknown  but  equal.  The  following  hypotheses  can  be  established: 


V  U1  =  ^2 

:  U1  4=  U2 

where  H0  =  null  hypothesis 

H-j^  =  alternate  hypothesis 

V>1  =  population  mean  vertical  velocity  of  severe  storms 
V*2  =  population  mean  vertical  velocity  of  non-severe  storms. 


The  formula  for  the  independent-t-test  is  (Hines  and  Montgomery,  1980) 

x  -  y  


t  = 


i  +  1 


n2 


Zx- 


-2 


-  nLx  +  Zy£  -  n2  y 


-2 


(100 


nl  +  n2  -  2 


with  n^  +  n2  -  2  degrees  of  freedom 

where  x^,  Yi  =  vertical  velocity  values  for  severe  and  non-severe 
cases  respectively 

nl»  n2  =  number  of  severe  and  non-severe  data  points 
x,  y  =  sample  means  for  severe  and  non-severe  cases. 


The  results  of  this  t-test  are  displayed  in  Table  XVI  when  using  data 


from  Figures  14  and  15. 


Table  XVI 


More  Statistics  on  the  24  April  and  6  May  1975  Storm  Complexes 
24  April  storm  Complex 

t0  =  2.63 
df  =  18 
a  =  .  05 

t. 025 ,  18  =  2.101 
tQ  >  t. 025 ,  18 

6  May  Storm  Complex 

tQ  =  2.2 
df  =  14 

using  a  95%  confidence  interval  (a  =  .05) 

t  n!  +  n2  -  2  =  t. 025 ,  14  =  2.145 
so  tQ  >  t.025,  14 

In  both  cases  above,  H0  is  rejected  and  there  is  a  95%  chance  that 
there  is  a  difference  between  severe  and  non-severe  thunderstorm 
updraft  velocities. 

Adler  and  Fenn  (1979  c)  examined  23  cases  of  storms  near  10  km 
from  both  24  April  and  6  May  1975.  Updrafts  from  severe  and  non- 
severe  thunderstorms  had  average  vertical  velocities  of  40.7  m  sec" 
and  19.92  m  sec“^  respectively.  Using  Equation  99,  severe  storms 
cooled  at  a  rate  of  2.4°k  min-^.  An  interesting  exception  to  these 
numbers  is  the  tornado  that  swept  through  heosho,  Missouri  at  0040 


Greenwich  time  on  25  April  1975.  Updraft  velocities  at  the  cloud  top 
reached  64.7  m  sec  ^  indicating  a  cooling  of  4. 0°k  min-^-  (Adler  and 
Fenn  1979  c) . 

Stereographic  Observations 

Hasler  (1981)  showed  that  stereographic  observations  of  cloud 
heights  and  their  temporal  changes  from  two  simultaneously  scanning 
visible  geosynchronous  satellities  is  a  fundamentally  new  tool  for  the 
atmospheric  scientist.  Stereographic  measurements  of  cloud  heights  are 
based  on  simple  geometric  relationships  and  boast  of  horizontal  reso¬ 
lution  of  0.5  km.  Conventional  measurements  with  geosynchronous 
satellites  have  a  spatial  and  temporal  resolution  of  1  km  and  3  min. 
Stereographic  capabilities  are  indeed  impressive  but  are  limited  to 
daylight  use.  A  shuttle-borne  lidar  only  requires  one  instrument 
and  can  be  used  in  day  or  night  with  horizontal  resolution  of  several 
meters . 

Mack  et  al  (1982)  analyzed  stereographic  data  from  several 
tropopause-penetrating  thunderstorms  over  Oklahoma  on  2  and  3  May 
1979.  He  determined  that  moderate  thunderstorms  grew  (cooled)  at  a 
rate  of  2°k  min  ^ .  These  results  agree  closely  with  that  of  Adler 
and  Fenn  (1979  c)  reported  in  the  previous  section. 

Damaging  Hailstorms 

A  large  number  of  damaging  hailstorm  producing  losses  near  $100 
million  struck  portions  of  Mew  Mexico,  Colorado,  Wyoming,  and  Montana 


from  May  through  August  1978.  Reynolds  (1980)  examined  these  storms 
with  visible  and  infrared  data  from  a  geostationary  satellite.  It  was 


determined  that  the  coldest  portion  of  the  storm  top  was  located  near 
the  area  of  hailfall.  This  cold  portion  of  cloud  top  was  1-8  degrees 
colder  than  the  surrounding  anvil  floor  which  was  located  at  the 
tropopause.  Reynolds  was  able  to  conclude  that  a  strong  correlation 
existed  "between  the  onset  of  large  hail  (>  25  mm)  and  cloud-top 
temperatures  becoming  colder  than  the  environmental  tropopause  tempera¬ 
ture"  (Reynolds,  1980:  345).  Table  XVII  depicts  typical  data  that 
Reynolds  worked  with. 

Reynolds'  discoveries  pertaining  to  hailstorms  implies  another 
use  for  a  shuttle-borne  lidar.  Instead  of  measuring  rates  of  change  in 
temperature,  the  lidar  could  take  several  measurements  across  the 
anvil  floor  within  the  80  sec  interval.  Those  towers  with  temperatures 
1-8°  coi<jGr  than  the  surrounding  top  would  be  identified  on  the  ground 
as  locations  of  possible  damaging  hail.  Figure  16  ideally  shows  a 
shuttle  lidar  making  two  measurements  of  a  storm  top  well  within  the 
80  second  limit. 

Scenario  II 

The  first  scenario  considered  in  the  previous  chapter  demonstrated 
that  the  shuttle  lidar  could  recover  a  single  temperature  measurement 
from  a  hypothetical  overshooting  top.  In  that  first  scenario,  the 
shuttle  was  orbiting  at  200  km  while  the  storm  top  was  at  10  km.  The 
retrieved  temperature  was  235°k.  This  second  scenario  will  continue 
the  first  and  measure  the  same  cold  area  at  a  later  time. 

It  is  now  known  that  the  lidar  can  track  the  same  target  area  for 
just  over  a  minute.  Assume  a  minute  after  the  first  measurement  the 


Table  XVII 

Damaging  Hailstorms  of  May-August  1978  (after  Reynolds 
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absorption  coefficient  UgQ  is  equal  to  .0438  km“l.  At  this  time,  the 

new  storm  top  temperature  and  height  are  unknown  except  that  the  top 

is  near  10  km.  The  value  of  Aq  from  Equation  41  will  still  be 
3/2  -l 

1742. 73°k  m  .  The  new  temperature  T*  can  be  found  still  through 
Equation  43  using  a  seed  guess  temperature  of  235°  (value  from  the 
first  measurement) 1 


T*i+1  = 


In 


* 


(1.439) (1421.436  cm-1) _ 

j'(1742.73°k3/2  m"1)  (.9999)  (235°k“3/2) 

(1  _  - 2 

'  (10) (.7486)  .0438  km 


=  230°k. 


So  in  one  minute,  an  energetic  updraft  drove  an  already  overshooting 
top  upwards  at  which  time  the  top  cooled  5°k  in  one  minute.  A  cooling 
rate  of  this  magnitude  would  classify  the  thunderstorm  as  intense 
(Mack  et_  al,  1982) .  The  storm  top  vertical  velocity  according  to 
Equation  96  is  .625  km  min~^. 

The  storm  top  is  now  located  10.625  km  above  the  surface.  The 
appropriate  military  and  civilian  organizations  would  be  notified  of 
the  severe  nature  of  the  above  storm  and  emergency  action  taken. 

The  80  second  loitering  time  per  storm  top  is  now  exeeded;  therefore, 


the  lidar  would  be  aimed  at  the  next  suspicious  thunderstorm. 


X.  Miscellaneous  Considerations 


Eye  Safety 

Table  V  (see  Chapter  VI)  shows  that  the  laser  considered  in  this 

thesis  develops  a  peak  power  of  1.8  MWatts.  This  is  a  great  deal  of 

power  and  would  damage  the  human  eye  if  the  laser  is  viewed  directly 

into  the  beam.  The  American  Conference  of  Government  Industrial 

Hygienists  in  1973  established  threshold  values  of  energy  density  for 

intrabeam  viewing.  A  pulsed  laser  operating  near  .7600  ym  with  pulse 

length  less  than  10*"^  seconds  has  a  threshold  energy  density  of 
-7  o 

6.5x10  J  cm""^  according  to  the  Conference  (Smith  and  Platt,  1977). 
Therefore,  there  would  be  no  eye  damage  if  laser  energy  density 
remains  below  this  threshold. 

The  space  shuttle  in  this  thesis  orbits  at  200  km.  It  must  be 
determined  if  the  eye  could  be  damaged  by  observing  the  shuttle  from 
the  ground  with  the  laser  operating.  The  shuttle  lidar  has  a  full 
angle  beam  divergence  of  1  mrad.  The  circular  area  illuminated  AL 
on  the  ground  can  be  found  by 

aL  =  W<ZL  tan  0/2>2  (101) 

where  ZL  =  height  of  shuttle 

6  =  full  angle  beam  divergence. 

Table  XVIII  clearly  shows  that  the  shuttle-borne  lidar  poses  no 
hazard  to  the  naked  eye  of  an  observer  on  the  ground.  This  will  not 
be  the  case  if  the  laser  is  viewed  through  a  high  powered  telescope; 
however,  this  shuttle  lidar  is  solely  used  to  remote  sense  temperature 


of  thunderstorm  tops  and  could  never  be  viewed  directly  from  the 


surface  due  to  cloud  attenuation. 


Table  XVIII 

Laser  Energy  Density  Reaching  the  Earth 1 s  Surface 


Pulse  Energy 
(J) 


Energy  Density  at  Surface 
(  J  cm~2) 


1 


3.18xl0~9 


10  3 . 18x10“® 

100  3 .18x10” 7 


Calibration  of  the  Lidar 

The  present  instrument  is  a  DIAL  system  featuring  overlapping 
co  linear  transmission  of  two  laser  beams.  DIAL  theory  according 
to  Equation  19  requires  that  the  extinction  properties  cf  the  atmos¬ 
phere  be  identified  with  exception  to  the  resonant  absorption  effects 
which  are  to  be  determined.  This  implies  that  the  instrument  must  be 
calibrated  so  that  both  lasers  react  in  a  similar  manner. 

Calibration  could  be  accomplished  by  tuning  the  on-line  laser  to 
the  off-line  frequency.  Simultaneous  measurements  from  both  lasers 
could  be  taken  of  transmitted  and  return  signals.  The  differences 
would  be  used  to  normalize  subsequent  data  runs.  The  on-line  laser 
would  then  be  tuned  back  to  the  absorbing  wavelength  and  operations 
would  begin.  Calibration  problems  have  been  eliminated  in  this  thesis 
by  assuming  atmospheric  extinction  at  both  wavelengths  i1-  the  same 
except  for  the  resonant  absorption  of  oxygen. 


Transmitter  Configuration 

Figure  17  is  a  block  diagram  of  a  lidar  transmitter  patterned  after 
the  DIAL  device  used  by  Kalshoven  et  al  (1981) .  Solid  lines  indicate 
the  passage  of  laser  light  while  dashed  lines  indicate  electrical  power 
lines.  The  krypton  gas  laser  is  used  to  pump  the  two  dye  lasers. 

Receiver  Configuration 

Figure  18  is  a  block  diagram  of  a  typical  lidar  receiver  used  for 
long  range  remote  sensing  taken  from  a  variety  of  sources  (Huf faker,  1978; 
Greco,  1980;  Kalshoven  et  al,  1981;  Spinhirne  et  al,  1982) .  This 
receiver  features  a  1.25  m  diameter  receiving  telescope  and  a  1  nm 
spectral  bandpass  filter. 

System  Arrangement 

Figure  19  portrays  an  idealized  lidar  system  arranged  in  the 
shuttle  cargo  bay.  Such  an  arrangement  can  easily  interface  with 
the  space  shuttle  or  Spacelab  (Greco,  1980) . 

Viewing  Geometry 

Figure  20  shows  possible  viewing  geometry  and  clearance  angles 
for  a  downward  viewing  lidar  (Huf faker,  1978).  When  considering 
the  clearance  angle  of  75  degrees  and  nadir  angle  of  60  degrees, 
the  lidar  could  search  a  maximum  area  of  3.8xl05  km2  while  operating 
from  a  200  km  orbit. 

The  Infrared  Equivalent 


This  thesis  has  attempted  to  show  that  a  laser  operating  from 
the  space  shuttle  can  measure  the  temperature  of  thunderstorm  cloud  tops. 


FLIGHT  ORIENTATION: 

T  PERPENDICULAR  ORBIT 
PLANE  (T^p) 

*1  NADIR 

-X  VELOCITY  VECTOR 


GROUND  TRACK 


At  this  point,  it  would  be  interesting  to  know  if  an  infrared  radio¬ 
meter  operating  from  the  shuttle  could  perform  the  same  feat. 

Adler  et  al  (1982)  flew  at  18  km  over  thunderstorms  in  NASA's  WB57F 
aircraft  equipped  with  a  11  urn  scanning  radiometer,  a  cloud  top  scanner 
(CTS).  The  CTS  was  accurate  to  approximately  +  2°K  at  temperatures 
colder  than  213°K,  The  instrument  also  features  an  80  m  resolution. 
Adler's  general  objective  was  to  compare  aircraft  derived  cloud  top  tem¬ 
peratures  witli  those  of  geosynchronous  satellites.  In  addition,  he  pro¬ 
vides  the  first  thunderstorm  cloud  top  thermal  mapping  indicating  the 
complex,  continually  changing  nature  of  bubbling  convective  cells.  It 
was  determined  that  for  typical  thunderstorms  the  satellite  derived  tem¬ 
peratures  were  about  15°K  warmer  than  the  aircraft  derived  temperatures. 
This  discrepancy  results  from  the  satellite's  large  field  of  view  (FOV), 

10x10  km,  so  that  the  temperature  determined  is  some  sort  of  average 

2 

value  for  the  100  km  area.  For  small,  growing  thunderstorms,  the 
satellite  over  estimated  top  temperatures  by  about  30-40°K  again  due  to 
FOV  effects.  According  to  Adler,  a  satellite  FOV  of  1  km  would  allevi¬ 
ate  these  problems. 

The  pioneering  work  of  Adler  et  al  (1982)  seems  to  indicate  that  it 
would  be  feasible  to  operate  a  CTS  with  a  FOV  of  1  km  and  measure  cloud 
top  temperatures  from  the  space  shuttle.  Nevertheless,  high  resolution 
infrared  scanners  generate  a  great  deal  of  data  which  requires  a  large 
processing  time  before  information  can  be  passed  onto  the  user.  Proces¬ 
sing  time  is  3  minutes  for  current  geosynchronous  infrared  scanners 
(Hasler,  1981).  The  shuttle-borne  lidar  is  not  a  scanning  instrument 
and  will  deliver  a  temperature  within  20  seconds  depending  on  the  sounding 


time  used. 


The  lidar  operating  at  200  km  achieves  a  horizontal  resolution  of 
10  m  for  a  storm  top  at  10  km.  It  is  doubtful  that  an  infrared  scan¬ 
ner  can  ever  be  designed  to  match  the  lidar* s  super,  high  resolution. 
The  shuttle  lidar  will  therefore  be  capable  of  conducting  a  finely  de¬ 
tailed  thermal  analysis  of  a  storm  top.  Such  an  analysis  could  indeed 
reveal  a  very  complex  cloud  top  temperature  structure  which  could  pro¬ 
vide  new  insight  into  thunderstorm  dynamics. 

The  pointing  of  an  orbiting  lidar  or  infrared  instrument  at  a  par¬ 
ticular  overshooting  top  is  a  minor  problem  that  could  be  solved  in 
several  ways.  First,  a  shuttle  qualified  payload  specialist  could  be 
trained  to  recognize  and  aim  the  instrument  at  a  promising  region  on 
the  cloud  top.  Second,  a  payload  specialist  could  coordinate  with  met¬ 
eorologist  on  the  ground  to  help  locate  a  suspicious  storm  top.  Ground 
based  personnel  would  use  radar  or  satellite  derived  data  to  aid  the 
payload  specialist  in  locating  a  top  that  deserves  special  attention. 


Ill 


XI.  Conclusion 


This  thesis  has  demonstrated  for  the  first  time  that  a  lidar 
operating  from  the  cargo  bay  of  the  space  shuttle  can  successfully 
measure  the  intensity  of  thunderstorms.  Thunderstorm  intensity  is 
monitored  by  measuring  the  time  rate  of  change  in  temperature  of  a 
thunderstorm’s  overshooting  top.  Severe  weather  occurring  on  the 
ground  takes  place  during,  or  just  after  a  period  of  rapid  cloud  top 
cooling. 

Temperature  is  measured  using  the  differential  absorption  lidar 
(DIAL)  technique.  This  technique  features  overlapping  colinear  trans¬ 
mission  of  two  laser  beams  at  different  wavelengths.  The  first  wave¬ 
length  is  set  at  .7696  urn  which  is  a  line  center  within  the  temperature 
sensitive,  absorbing  portion  of  the  oxygen  "A"  band  centered  near  .7700 
urn.  The  second  laser  wavelength  is  set  at  .7614  um  which  is  a  nearby 
spectral  position  where  absorption  from  the  oxygen  HA"  band  is  negli¬ 
gible.  This  DIAL  technique  is  designed  so  that  scattering  and  absorp¬ 
tion  properties  of  the  atmosphere  are  identical  except  for  oxygen  res¬ 
onant  absorption  at  line  center.  It  is  shown  that  the  lidar  can  measure 
the  oxygen  resonant  absorption  coefficient  at  the  heights  of  typical 
thunderstorms.  A  temperature  value  is  recovered  by  placing  the  retrieved 
absorption  coefficient  into  a  quickly  converging  iterative  expression. 

The  two  lasers  required  by  the  DIAL  technique  are  dye  lasers  which 
are  optically  pumped  by  a  krypton  gas  laser.  The  laser  system  output 
is  a  conservative  300  watts  which  is  within  space  shuttle  power  con¬ 


straints 


A  variety  of  hypothetical  thunderstorms  are  probed  ranging  in 
height  from  10-17  km  at  different  shuttle  altitudes  ranging  from 
100-250  km.  Signal-to-noise  calculations  are  performed  considering 
the  worse  case  at  day  with  the  sun  directly  overhead  and  at  night 
with  a  full  moon  directly  overhead.  Success  of  the  system  is  based 
on  useful  values  of  the  signal-to-noise  ratio. 

The  shuttle- borne  lidar  will  remain  within  range  of  a  given 
overshooting  top  for  80-100seconds  while  orbiting  at  200  km.  During 
this  time,  the  rate  of  change  in  temperature  can  be  determined  by 
taking  two  measurements  of  the  same  overshooting  top  within  a  minute 
of  each  other.  A  storm  is  considered  intense  or  severe  if  the  rate 
of  change  in  temperature  is  greater  than  2°K  min-^.  Alternately, 
during  the  80-100  second  loitering  time,  4  single  temperature  measure¬ 
ments  could  be  made  of  separate  overshooting  tops.  Those  tops  determined 
to  be  1-8  degrees  colder  than  the  tropopause  correspond  to  locations 
on  the  ground  where  damaging  hail  (£25  mm  in  diameter)  can  be  expected. 

A  shuttle- borne  temperature  sensing  lidar  would  be  of  great  re¬ 
search  value  to  the  atmospheric  sciences.  Horizontal  resolution  is 
10  m  for  a  storm  top  at  10  km  and  shuttle  operating  at  200  km.  This 
implies  that,  for  the  first  time,  high  resolution  thermal  mapping  of 
storm  tops  can  be  conducted  from  orbit,  A  shuttle-borne  lidar  could 
supplement  existing  severe  storm  monitoring  networks  by  providing  timely 
warnings  of  intensifying  storms.  In  this  way,  such  a  sensor  could  pre¬ 
vent  millions  of  dollars  in  property  damage  and  possibly  save  lives. 


Re  c  ominend  a  t  i  on  s 


As  an  effort  to  guide  future  research  in  the  area  of  shuttle- 
borne  lidars  used  to  measure  the  intensity  of  thunderstorms,  some 
recommendations  can  be  made.  The  most  obvious  improvement  that  could 
be  made  is  to  investigate  increasing  the  power  available  to  the  lidar. 
Nuclear  power  would  not  be  needed.  Additional  fuel  cells,  each  weighing 
747  kg  and  delivering  840  kWh  (Huffaker,  1978),  could  be  placed  in  the 
shuttle's  cargo  bay. 

Investigations  could  be  conducted  into  increasing  the  pulse  rep¬ 
etition  rate  (PRR)  of  dye  lasers  operating  in  the  near  infrared  portion 
of  the  spectrum.  As  the  PRR  increases,  small  increases  in  the  signal- 
to-noise  ratio  can  be  realized. 

Chapter  V  mentions  that  a  cirrus  canopy  resulting  from  previous 
convection  will  sometimes  form  above  a  thunderstorm  thus  masking  the 
storm  top  from  an  orbiting  sensor.  This  situation  was  not  considered 
here  but  future  work  could  attempt  to  model  laser  beam  penetration  of 
the  cirrus  canopy  before  the  beam  intercepts  the  storm  top  to  determine 
if  a  useful  signal  can  be  recovered.  Perhaps  actual  experimentation 
will  be  required. 

Finally,  current  theory  claims  that  a  lidar  can  measure  the  at¬ 
mospheric  state  variables  (wind,  humidity,  pressure,  and  temperature 
according  to  Atlas  and  Korb,  (3981)).  It  would  be  interesting  to 
investigate  developing  an  economical  shuttle-1 idar  system  which  would 
be  capable  of  measuring  the  atmospheric  state  variables.  This  in¬ 
novation  would  lead  to  a  complete  orbiting  weather  station. 
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